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PREFACE 

The data and data formats described herein have been generated in cooperation 
with several people at the Lockheed Martin Advanced Technology Center, formerly 
the Lockheed Research and Development Division, in particular with R. D. Sharp, 
E. G. Shelley, W. K. Peterson, D. L. Carr, and W. E. Francis. The archiving of 
these data was supported by NASA Contract NAS5-33047 and Grant NAGW-4496. 

The magnetic fields included here are based on ISEE-1 magnetometer data made 
available by C . T. Russell of the Institute for Geophysics and Planetary Physics 
at the University of California, Los Angeles, California. The processing and 
editing of the magnetometer data are the sole responsibility of the author. 

Any questions regarding these data files should be directed to: 

Dr. 0. W. Lennartsson 
Lockheed Martin Missiles & Space 
Advanced Technology Center 
Dept. Hl-11, Bldg. 255 
3251 Hanover St. 

Palo Alto, CA 94304 

Telephone: (650) 424-3259 

Telefax: (650) 424-3333 

e-mail: lenn0space.lockheed.com 
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1. INTRODUCTION 


The ISEE-1 spacecraft {along with the ISEE-2) was launched on October 22, 1977, 

into an orbit with apogee at almost 23 RE (geocentric) , perigee at about 300 km 
altitude, an inclination of 29 deg, and an orbital period of 57 hours. It was 
placed in a spinning mode with the axis nearly perpendicular to the solar eclip- 
tic plane and with a period of approximately 3 seconds. It reentered almost 10 
years later, on September 26, 1987. 

The Lockheed Plasma Composition Experiment flown on ISEE-1 is one of a family 
of instruments using the same type of ion optics and covering nearly the same 
range of energies (0 eV/e to about 17 keV/e) which have also been flown on GEOS 
1 and 2, Dynamics Explorer One, and AMPTE CCE [Shelley et al., IEEE Trans. 
Geosci. Electron., GE-16, 266, 1978]. The Principal Investigator of the ISEE-1 
experiment during its development and during its several years of operation in 
orbit was R. D. Sharp. In the later phase of data analysis this role was trans- 
ferred to O. W. Lennartsson. For a review of the first few years of data anal- 
ysis see Sharp et al. in Energetic Ion Composition in the Earth's Magnetosphere, 
edited by R. G. Johnson, TERRAPUB, Tokyo, 231, 1983. 

The ISEE-1 instrument consists of two nearly identical mass spectrometers with 
the respective fields of view centered 5 deg above and 5 deg below the spin 
plane, that is about 5 deg above and below the solar ecliptic plane. Each field 
of view is 10 deg wide along the spin plane, and 10 to 50 deg wide transverse to 
this plane, being the widest at the low energy end and gradually decreasing to 
10 deg with increasing energy. All data used for this archiving were obtained 
with one spectrometer, the one viewing below the spin plane. Information on the 
instantaneous pitch angles is provided by the ISEE-1 Fluxgate Magnetometer 
[Russell, IEEE Trans. Geosci. Electron., GE-16, 239, 1978], 

Each spectrometer consists of an electrostatic analyzer to select energy per 
charge, followed by a combined electrostatic and magnetic analyzer to select 
mass per charge. Both analyzer sections have particle detectors, so at each 
energy setting the experiment provides both the total ion flux and the partial 
flux at a selected mass per charge. The detector in the energy analyzer is off- 
set from the center path of ions and corresponds to a view direction that is 
slightly different from the nominal mass spectrometer field of view, being more 
nearly within the spin plane. To be more specific, the two detectors have about 
the same size fields of view at all energies, but the "energy detector" central 
view direction varies, being some 15 deg away from the spin plane on the oppo- 
site side of the "mass detector" view direction at near-zero energies, and grad- 
ually approaching the spin plane with increasing energy. 

Each combination of energy and mass is maintained for at least 1/16 sec in high 
telemetry bit rate and 1/4 sec in low (normal) bit rate. Different combinations 
are stepped through in a cyclic fashion according to various patterns controlled 
by a random access memory which is programmable from the ground. Some patterns, 
or "modes", may require only a couple of minutes per cycle, others may take 15 
to 20 minutes. 

The mass selections usually include one that blocks all ions from reaching the 
second detector, which is a modified Johnston electron multiplier, allowing 
intermittent measurements of the noise associated with penetrating radiation 
(mostly MeV electrons and associated bremss t rahlung) . These measurements are 
later used to correct the count rates of mass analyzed ions. The background 
count rate in the first detector, which is a so called spiraltron, is generally 
negligible and is not corrected for. 



The maximum energy range is 0 eV/e (or spacecraft potential) to 17.9 keV/e, 
divided into 32 contiguous channels, but the lowest channel, from 0 eV/e to 
about 100 eV/e, is normally limited to energies above 10 eV/e by an RPA (re- 
tarding potential analyzer) in the entrance. That same RPA is used to provide 
"cold plasma" data from 0 to 100 eV/e (retarding within the lowest channel) 
during part of some measurement cycles. Because of measurement uncertainties 
associated with spacecraft charging and plasma convection, the "cold plasma" 
data have not been included here. The lowest "hot plasma" channel, the one 
covering 10 eV/e to 100 eV/e, has been used only sparingly for those same rea- 
sons . 

All velocity moments, with one exception (see below), use 100 eV/e as the lower 
energy limit. As far as moments are concerned, the highest energy channel 
(above 16 keV/e) has not been used either, because of a slight variation over 
time of the mass peak locations in that channel. The energy range of moments 
is therefore generally 100 eV/e to (about) 16 keV/e, unless further limited by 
a particular instrument energy/mass scan mode (see below) . 

The mass range is from below 1 amu/e to about 150 amu/e, divided into 64 chan- 
nels, but the most common scan cycles only sample 4 to 5 channels, corresponding 
to the principal ions, plus a background channel (below 1 amu/e) . 


2. ARCHIVAL DATA FORMAT 

These archival da£a are a combination of ion data, magnetic field data, and s/c 
ephemeris, generally sorted by instrument energy/mass scan cycle and placed in 
five kinds of ASCII files, called the ED-, MD-, MC-, MS-, and EX-files, respec- 
tively. There is one set of these files for each day of data, and the date is 
part of the file names, represented by the last two digits of the year, followed 
by three digits for the day of the year, as for example ED78057.DAT, etc. The 
format and content of each kind of file are described in detail in following 
sections . 

The ion data are a combination of velocity moments, calculated two different 
ways (using two different sets of assumptions about symmetry in velocity space) , 
and energy and mass spectra. The velocity moments are calculated both from the 
total ion count rates in the energy analyzer, assuming that all counts are due 
to H+ ions (and using one set of symmetry assumptions) , and from the partial 
count rates at M/Q - 1, 2, 4, and 16 in the mass analyzer, assuming those counts 
to be due to the four principal magnetospheric ions H+, He++, He+, and 0+, re- 
spectively (and using both sets of symmetry assumptions) . The "total ion" mo- 
ments are listed in the ED-file, the partial (mass analyzed) moments are listed 
in the MD-file. 

The energy spectra are based on mass analyzed counts only, and are listed in the 
MC-file for the same four principal ions (same four M/Q values) , plus a fifth 
ion, namely the 0++ (M/Q= 8) . The five energy spectra are a compacted represen- 
tation of the raw count rates at each energy (spin-averaged and maximum) . 

The mass spectra (when produced) show the accumulated raw counts (and number of 
samples) in each of the 64 mass channels in four energy ranges. These are list- 
ed in the MS-file. 

The magnetic field data (time averaged) are listed in all of the ED-, MD-, MC-, 
and MS-files. In addition, the magnetic field measurements have been used in 
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conjunction with the ion data to determine pitch angles and ion plasma beta. 


The s/c ephemeris is also listed with the ion data in the ED-, MD-, MC-, and MS 
files, but an extra set of half-hour interval parameters are listed in the EX- 
file. That file contains only ephemeris, and has been used to label certain 
data plots. These plots are available in hardcopy (see end of this document fo 
description) . 


3. TIME COVERAGE 

These data have been produced to cover the entire 4.4-year lifetime of the 

Lockheed ISEE-1 instrument, from launch on October 22, 1977, through March 20 

of 1982. The coverage is continuous with the following exceptions: 

a. The initial two weeks were devoted to extensive in-flight testing and 
calibration of the experiment. The resulting data have not been used 
for this archiving, so the beginning date is day 312 of 1977. 

b. The instrument was normally turned off for 1.5 to 2 hours at perigee, 
occasionally longer (primarily to avoid the inner radiation belts) . 

c. Extended periods of "cold plasma" measurements (a few tens of hours 
all together) . No such measurements are included here, but they were 
normally carried out as brief (few minutes) interruptions of longer 
energy /mass scan cycles and do not normally show up as data gaps . 

* 

d. Solar wind measurements in the so called "sun-synchronous" mode. 

In this mode data were taken at an increased rate in a narrow spin 
angle sector around the sun direction, and no data were taken at 
other angles, in order to maintain the same average data rate per 
spin. These data have been judged too cumbersome to treat within 
the particular format of this archiving. The resulting data gaps 
are numerous and extensive, sometimes leaving an entire day blank. 

e. Gaps in the data transmission/reception or poor quality reception. 
These are usually of modest duration (few hours) but fairly common. 


4. NUMERICAL ACCURACY 

The absolute calibration of the instrument is believed to have been accurate to 
better than 30% in the early phase of in-flight operations. This estimate is 
based mostly on pre-flight tests in the laboratory, but also on extensive early 
in-flight confirmation of the location and shape of major peaks in M/Q spectra. 
As discussed below (in connection with the MC-file) , the first 28 months of 
mass analyzed data have been derived from counts acquired at the higher of 
two detector triggering levels, in order to minimize the effect of background. 
This triggering level has become progressively less sensitive over the years, 
but it has been monitored, and the count rates, through February (actually 
through March 1) of 1980, have been adjusted accordingly, before being used in 
the moment calculations. These adjustments are believed to have maintained the 
bsolute calibration accuracy to within about 30% to 35% during high (normal) 
sensitivity operations. Data taken after February of 1980 are based on the 
lower triggering level, which is believed to have remained fairly constant over 
time. Except for the somewhat higher background in these later data, the abso- 
lute accuracy is expected to be comparable to that of the earlier data. 



4.1 Normal High Sensitivity Mode 


These estimates assume ideal conditions, however, which means that the differen- 
tial flux of any ion species is isotropic within the field of view of each sam- 
pling, constant in energy over the width of each energy channel, and constant in 
time during each sampling. Any or all of these conditions may be violated in a 
real situation, and the resulting error may be significantly greater than 35%. 

To take an extreme example: if a phase space density is inferred from the count 

rate in the lowest energy channel, from between 10 eV/e and slightly beyond 100 
eV/e, that count rate is divided by 40 eV/e (the weighted center energy) . If in 
fact the ions all have energies near the lower or upper edges of this channel, 
then the phase space density will be either four times too small or almost three 
times too large. Errors of a similar magnitude may occur when the ion flux has 
an extremely narrow angular distribution, or when it varies rapidly in time. 

The energy distribution is of less concern at higher energies (smaller delta-E 
over E) , but the potential error in phase space density due to narrow angular 
distributions can be large at any energy, especially in low (normal) bit rate, 
where the instrument sweeps through 30 deg of spin angle during a sampling 
(1/4 sec) . 

Summing (weighted) phase space densities into moments does reduce the relative 
measuring errors, but to infer complete three-dimensional moments from these 
data requires assumptions about that part of the unit sphere which is not being 
sampled. The potential new errors are impossible to calculate, of course, but 
they can be minimized by making realistic assumptions. The moments listed in 
the MD-file have been calculated two ways, using two entirely different sets of 
assumptions (see file description) , each set being optimal in some part of the 
ISEE-1 environment. As can be seen from the listed numbers, the two ways of 
calculating number densities, for example, produce nearly the same result much 
of the time, at least in the magnetosphere. If the respective set of assump- 
tions is correct, then the moment calculations are intrinsically accurate to 
within a few percent. 


4.2 Low Sensitivity Mode 

If operated normally, the particle counters will saturate when exposed to the 
high directional flux of H+ in the magnetosheath and solar wind. Even the He++ 
may cause saturation in these regions of space. To avoid this and to increase 
the dynamical range of the data, the instrument was designed to have a second 
commandable detector function that would reduce the sensitivity by about two 
orders of magnitude. This is accomplished in one case by turning off one of two 
power supplies (in the mass analyzer) and by shifting the power to another de- 
tector in the other case (in the energy analyzer) . For the most part, this sen- 
sitivity change is programmed to occur automatically when the count rate exceeds 
a certain threshold, and to be deactivated at the start of the next group of 
commands, that is typically within 8 sec in low bit rate and within 2 sec in 
high. New commands are activated at intervals that are multiples of 32 sam- 
plings, called "octaves", usually at every such octave. Normally the instrument 
maintains the same mass channel for one or more whole octaves at a time as well. 

In theory this dual sensitivity function is simple enough, but it has several 
practical problems, especially when it is set to be automatic. The change to 
low sensitivity takes a finite known time (close to one second) , but the telem- 
etry bit stream does not allow the exact starting time to be shown, only the af- 



fected octave is flagged. Instead the time of change must be estimated from the 
sequence of count rates, which is hampered by rapid strong variations in the ion 
flux. In automatic operation it is also reset to normal in an automatic fashion 
(takes negligible time), so the change to low sensitivity may be repeated as of- 
ten as every octave (or every other octave in high bit rate) . 

During the computer processing of these archival data, the transitions from high 
to low sensitivity have been treated as data dropouts, for simplicity, and only 
counts taken at least 1.5 seconds (1/2 spin) beyond the estimated start of each 
transition have been included, after multiplication by a constant scaling factor 
representing the end state (no such multiplication in the MS-file) . This method 
avoids having the transient count rates multiplied by the wrong variable scaling 
factor, due to uncertainties in the exact timing, but it creates a systematic 
undersampling of the ion flux near its maximum and on one side of this maximum. 
The situation is further aggravated by the fact that these automatic sensitivity 
changes tend to occur where the ion fluxes have the strongest anisotropies and 
the narrowest energy distributions. 

In view of these difficulties, all data obtained in low sensitivity mode should 
be treated with suspicion. They show the gross and qualitative plasma proper- 
ties very well, but the tabulated numbers may be substantially wrong. If it is 
important to establish absolute values, it is worth comparing these numbers with 
archived data from some of the other particle instruments on the ISEE-1 and -2. 
The data most often obtained in this mode are those in the ED-file and the H+ 
data in the MD-, MC-, and MS-files. Sometimes the He++ data are also affected. 
All data obtained this way are clearly flagged (see file descriptions) . 


5. KNOWN ERRORS 

All data, including the magnetic field data, have been screened very rigorously 
with respect to potential problems in the telemetry transmission and reception. 
Any time the telemetry quality flags suggest a potential problem, all data that 
could possibly be affected have been discarded without further investigation. 

In addition, every ion count rate and every magnetic field reading have been 
examined for clearly "unrealistic" values, that is unrealistically large values 
under given circumstances (given region of space) . Such values have been dis- 
carded on an individual basis (sample by sample) . 


5.1 Ion Data 

The remaining errors, as far as they are known, are rather few and normally 
obvious from the wider context, especially when the various files are intercom- 
pared. One group of known errors consists of a few isolated "bursts of counts” 
in the MS-file (mass spectra) at M/Q values where no significant ions are to be 
expected. Being outside of the nominal mass peaks, these false counts are not 
included in the MC- and MD-files. Another group of suspected errors consists 
of a few isolated "bursts of ions" having strangely repetitive count rate and 
angle of motion, as evident in the MC-file. 

Another known but subtle error may occur in the lowest energy channel (< 100 
eV/e) in mass analyzed H+ data, at very low count rates. This error is most 
likely caused by scattered sunlight. Specifically, when the instrument aperture 
is pointing very nearly sunward, the energy analyzer is set for E/Q < 100 eV/e, 
and the mass analyzer is set for M/Q = 1, a very weak signal may arise from what 
is believed to be ions produced internal to the instrument preacceleration re- 



gion by photoionization of residual gas, giving the impression of a weak anti- 
sunward flow of H+ ions. No systematic investigation has yet been undertaken 
to quantify this problem, but the effect is probably only of concern in the tail 
lobes, where natural H+ counts may be extremely low [Sharp et al . , J. Geophys . 
Res., 86, p. 4639, 1981]. This is one of the reasons why the lowest energy 
channel is treated separately in the moment calculations (see below) . 


5.2 Magnetic Field Data 

These have more errors in them, unfortunately, but that is offset by the fact 
that expertly processed values can be obtained from other archival files sup- 
plied by the Principal Investigator for the ISEE-1 Fluxgate Magnetometer (C. T. 
Russell; NSSDC identification No. 77-102A-04; see also newly released data on 
the World Wide Web, at http://www-ssc.igpp.ucla.edu/forms/isee/). There are 
two main types of errors occurring here, one brief error often present at the 
gain changes (at most once or twice per orbit) due to a misalignment of the gain 
flag and the actual data on the tape, and another less common but longer lasting 
error that may appear in connection with certain offset corrections (related to 
the antenna flip status) . These errors are more fully explained in the descrip- 
tions of the ED- and MD-files, along with methods to recognize the errors (see 
also description of hardcopy plots at the end of this document) . Errors in the 
magnetic field only affect part of the ion data, namely the beta value and the 
moments that are based on pitch-angle binning (the second method used in the MD- 
file) . 


6. THE ED-FILE 

This file contains data from a particle counter called the "energy detector", 
or "ED", which intersects the ions after energy analysis but before mass analy- 
sis. The count rates have been converted to certain velocity moments assuming 
that all counts are due to protons . These moments are thus not directly compa- 
rable to any of the moments in the MD-files, unless the protons are indeed the 
dominant species. Even in that case, there may often be significant differences 
between the ED and MD proton moments, because the ED and MD particle counters 
respond to ions with slightly different external angles of motion (see above) , 
and the two sets of count rates are usually averaged over different times. 

Since the ED count rates are independent of the mass channel selection, and the 
energy channels are normally scanned at a higher rate than are the mass chan- 
nels, it is usually possible to obtain several consequtive sets of ED moments 
during each complete instrument scan cycle. This advantage has been utilized 
whenever possible, and the ED averaging intervals are more or less independent 
of the phase of the instrument cycle. The criterion most commonly used for de- 
fining the ED intervals is to have a "sufficiently dense" coverage in energy 
and spin phase angle. This criterion varies with different types of energy and 
mass scan patterns, but typically amounts to having at least every other energy 
channel sampled in each of 12 spin angle sectors (30 deg wide sectors) . The 
missing matrix points are filled in by interpolation between adjacent points 
before the moments are calculated. 

The criterion for ending each ED averaging interval is modified when the instru- 
ment is in low sensitivity mode (usually set by automatic triggering at extreme 
count rates) and, at the same time, the count rates show strong (10-fold or 
stronger) spin phase modulation. At such times the ED intervals are made equal 
to the instrument scan cycle. This situation occurs most often in the solar 



wind. The purpose is to suppress artificial oscillations in the moments caused 
(in some modes) by spin phase dependent energy samplings. The modification has 
not been applied uniformly throughout the entire data set, however, so there are 
periods of "strange" and non-physical modulations in the ED moments, especially 
in the densities. At any rate, all data taken in low sensitivity mode should 
be treated with care; they are not intended to represent accurate quantitative 
measurements, but they have been included to complete the gross qualitative 
picture and to provide continuity between times of normal magnetospher ic data. 
All low sensitivity data are flagged as such (see below) . 

Unlike the MD data, the ED data have not been corrected for background counts 
due to penetrating radiation (not measured) , but they are known to be much less 
affected by that problem than are the MD data (different kinds of detectors) . 

The ED data also differ from the MD data in that no standard deviations are cal- 
culated from the combination of counts (a saving in CPU time and output volume) . 
If necessary, approximate upper limits on the ED standard deviations can be es- 
timated from the corresponding standard deviations in the MD-file (the MD data 
have additional variance associated with background correction) . 


6.1 File Format 

The ED-file has been written by a formatted sequential FORTRAN WRITE: 

WRITE (9, 210) IYYDDD, IUTSEC, BXED, BYED, BZED, BTED, ITHETB, IPHIB, 

*EDDNS, EDEMN, EDVEL, IEDANG, ACSEC , XCSEC , ENRG1 , ENRG2 , RX, RY, R2 

210 FORMAT (216, 4F8 .1,14,15,3 (1PE10 .2) ,15,2 (0PF8 .2) , 2F5 . 1, 3F6. 1) 

The variable names represent the following quantities: 

IYYDDD is the year (two digits) and day of year (three digits) . 

IUTSEC is the universal time in seconds at the midpoint of the averaging 
interval. The nominal length of each interval is defined by the 
distance between adjacent midpoints, but the actual length may vary 
somewhat because of time gaps (noise) in the data (see ACSEC below) . 

BXED, BYED, BZED, and BTED are the GSE components and absolute value, respec- 
tively, of the measured magnetic field, each averaged over the ED 
interval (for data quality check, see next section) . The unit is 
"nanotesla" ("nT"), or equivalently, "gamma". 

ITHETB and IPHIB are, respectively, the magnetic field elevation and longitude 
angles in GSE coordinates. ITHETB ranges from -90 deg (southward) to 
+90 deg (northward) , IPHIB from -180 deg (antisunward) to +180 deg (anti- 
sunward) , with 0 deg being sunward and 90 deg duskward. These angles 
have been calculated from the time averaged field components, and are 
rounded to the nearest integer. 

EDDNS is the ion number density in "/cm3" (assuming protons) contained in the 

nominal energy range from 100 eV to about 16 keV (see also next section) . 

EDEMN is the "thermal" ion energy in "keV", that is the mean energy (0.1 - 
16 keV) minus the energy associated with the common drift motion. 

EDVEL is the common drift speed in "km /sec" in a plane parallel to the s/c 



spin plane, that is approximately in the solar ecliptic plane. 

IEDANG is the longitudinal direction angle of that same drift, with 0 deg 
being sunward, and 90 deg duskward (-180 deg to +180 deg) . 

ACSEC is the number of seconds of data included in the moments. This number 
is normally at least 6% smaller than the length of the averaging inter- 
val, since neither the lowest nor the highest energy channels are in- 
cluded (may be further reduced by data gaps) . 

XCSEC is the number of seconds during the averaging interval when the ED detec- 
tor was in low sensitivity mode. This number is almost always zero in 
the magnetosphere, but usually nonzero in the magnetosheath and solar 
wind. If nonzero, it implies that at least some of the count rates have 
been multiplied by a crude scaling factor {= 130) to compensate for the 
reduced sensitivity. 

ENRG1 and ENRG2 are, respectively, the lower and upper energy limits used in the 
moment calculations, in units of M keV" . These are normally 0.1 and 16.1, 
but may be modified in some modes (cf . comments on energies above) . 

RX, RY, and RZ are the GSM coordinates of the s/c (at midpoint of interval) in 
units of "earth radii”. On the rare occasions when the GSM coordinates 
are unavailable (due to faulty ephemeris tape) the GSE coordinates are 
used instead (available on raw data tape) . Those occasions are flagged 
in the EX-file (also flagged in the MC-, MD- and MS-files) . 


6.2 Further Explanations of Certain Variables. 

The magnetic field values provided here may to some extent duplicate what is 
now available at the NSSDC (Prof. C. T. Russell; NSSDC identification number 
77-102A-04) and on the World Wide Web (http: //www-ssc .igpp.ucla . edu/ f orms/isee/ ) 
but are not intended to supplant those. The values derived here may not have 
been adequately screened against telemetry noise or adjusted for instrument 
anomalies, so it is recommended that the PI provided data be consulted whenever 
there is doubt about accuracy. There are two types of errors known to occur 
here : 

a. When the magnetometer is commanded to low gain on the inbound leg of 
the s/c orbit, usually between R = 8 Re and ~ 5 Re, the corresponding 
flag may show up slightly too late on the raw data tape, resulting in 
a brief underrepresentation of the field by a factor of 32. This may 
last for about a minute or less (part of one major frame) and affects 
a single ED interval. This error ought to be fairly obvious in the ED 
data, and therefore traceable in the MD data as well (by comparing the 
times) . The corresponding mismatch at high gain command on the out- 
bound leg, that is a sudden 32-fold increase of the field, is usually 
discovered as "unreal” and corrected. 

b. Part of the computing procedure to decipher the raw telemetry data is to 
infer, indirectly, the current magnetometer antenna flip status from the 
uncorrected field values themselves (pertains to the s/c Y and Z coordi- 
nates) . This part fails occasionally, which in turn makes the remaining 
procedure unable to correct the field readings for magnetometer offset, 
resulting in anomalous modulations of the GSE components. Early in the 
mission, when the offset is fairly small, these modulations cause the 



sum BXED**2+BYED**2+BZED**2 to be much smaller than the square of the 
time-averaged absolute field strength BTED**2. Later in the mission,, 
when magnetometer drift has created a larger offset, this sum may settle 
at a rather large artificial value, of order 100 nT. In either case, it 
is recommended that the two measures of field strength be compared rou- 
tinely, and that other magnetic field data be consulted if those two 
differ by more than a few percent. This problem usually only occurs in 
regions of weak magnetic field, mostly in the solar wind, where it is 
easily recognized by the aforementioned sum being nearly constant, with 
no evidence of diamagnetic influence by the ion pressure (see below) . 

An additional (or alternative) check can be made using certain MD data 
(see below) . 

The velocity moments have been calculated from a 32 energy channel (covering the 
entire range from 10 eV to about 18 keV) by 12 spin angle sector (30 deg each) 
matrix of time-averaged count rates in the following steps: 

I. A phase space density is assigned to each matrix point, using the 
local count rate when available, or interpolating between adjacent 
count rates if the point has not been sampled. The bottom and top 
energy channels are included in the interpolation procedure, when 
necessary, but not in the integration (summation) over energy. 

II. Within each energy channel the phase space densities are weighted 
by cosines or sines of the spin angle (center of sector) and summed 
over angle to form two orthogonal projections. These projections 
are in turQ weighted by energy and by an energy bandwidth (see below) 
and summed over energy channels, forming two orthogonal components, 
approximately the GSE X and Y components, of a vector that is pro- 
portional to number flow density. Both components, are then divided 
by a total (scalar) sum of phase space densities weighted by the 
energy bandwidth and by the square root of the energy, that is by a 
sum proportional to number density, to form (approximate) X and Y 
drift velocity components, and an angle IEDANG . 

III. A spherical coordinate system is envisioned with its polar axis along 
the drift velocity vector, that is in the GSE X-Y plane. It is now 
assumed that the phase space density has azimuthal symmetry in this 
coordinate system (rotational symmetry around drift velocity vector) . 

The number density, drift speed (but not angle), and mean energy are 
recalculated by summing over solid angle and energy in this system. 

The solid angle weighting factors in this case are zones on a unit 
sphere, each defined at the intersection with the GSE X-Y plane by the 
boundaries of a 30 deg spin angle sector, or by one boundary and the 
drift velocity vector. These zones are partially overlapping, and 
summing over 12 spin angle sectors (typically * 14 zones) means cover- 
ing the unit sphere twice, so a factor 1/2 is applied to each sum. 

The drift speed (EDVEL) and mean energy (EDEMN) are obtained by divid- 
ing flow density and energy density, respectively, by the number den- 
sity (EDDNS) . The mean energy is converted to "thermal" energy by 
subtracting the equivalent drift energy. 

If the drift speed in step II is less than 14 km/sec (less than 1 eV energy), 
the coordinate system is instead aligned with the X-Y projection of the magnetic 
field, that is with the angle IPHIB, provided the elevation angle ITHETB is be- 
tween -45 and +45 deg. If the latter is not the case, and if the drift speed is 
below 14 km/sec, the coordinate system is aligned with the s/c spin axis, and 



the phase space densities are treated as isotropic. Only density and mean en- 
ergy are recalculated in these two cases. 


The summation over energy treats each energy channel, except the second one 
(first channel not included), as a point measurement at the center energy, and 
takes the energy bandwidth to be the distance between adjacent channels in the 
trapezoidal fashion. At the second channel an extra term is added to extend the 
integral from the center energy downward to about 0.1 keV, assuming the flux to 
be a constant. This addition brings the mathematical energy range in better 
agreement with the instrumental range of acceptance. 

The variable ACSEC is a sum of elementary time segments associated with each 
commanded setting of the power supplies controlling energy and mass channels. 
These elementary times are 1/4 sec during low bit rate operation (about 80% of 
the time) and 1/16 sec during high bit rate operation. These times include the 
resetting of the power supplies, however, and are slightly longer than the times 
associated with particle counting. For simplicity, the particle counting is 
interrupted for about 12% of the elementary time segments in both low and high 
bit rate operation to allow for resetting. The ACSEC therefore exceeds the ac- 
tual particle counting time by about 14% in both cases. 


7. THE MD-FILE 

This file contains data from a particle counter called the "mass detector”, or 
"MD", which receives the ions after both energy and mass analysis. Count rates 
of the four principal magnetospheric ion species, the H+, He++, He+, and 0+, 
have been sorted by 32 energy channels (10 eV/e to 18 keV/e), 12 spin angle sec- 
tors (30 deg each) and 9 pitch-angle ranges (20 deg each) and averaged over each 
complete energy /mass scan cycle of the instrument. At the end of each cycle the 
averaged count rates have been converted to certain velocity moments using two 
different assumptions: the velocity distributions have rotational symmetry a- 

round either (A) the bulk flow vector or (B) the local magnetic field direction 
(see below). Note: the ion labels are applied to certain M/Q values (the in- 
strument does not measure mass per se) and may be inappropriate in the magneto- 
sheath and solar wind, especially at M/Q= 4 (He+) . 

The energy/mass scan cycles vary in length from a few minutes to about 20 min 
or more, depending on the instrument mode of operation. Each cycle may cover 
only a few mass channels and a reduced energy range, or it may cover the full 
energy range at each of 64 mass channels (load mode) . The most common modes in 
the magnetosphere provide for multiple energy scans at each of 5 to 7 mass chan- 
nels (including a background channel at M/Q < 1), where each energy scan may 
sample a different subset of the 32 energy channels (e.g. every fourth channel 
in four interleaving scans) , and each energy channel is maintained for about 
four seconds (1.3 s/c spin cycles). In addition, each cycle may contain a few 
scans through all 64 mass channels at a few energies, as well as one or two 
brief scans through the RPA voltages in the lowest energy channel (cold plasma 
measurements not included here) . In order to simplify tabulation of moments, 
no distinction is made here between different phases of a given cycle. That is, 
all moments are treated as averages over the same cycle, although different ions 
were in fact sampled at different phases (the same ion often more than once per 
cycle) . 

All MD moments have been corrected for background counts due to penetrating 
radiation (mostly MeV electrons and associated bremsstrahlung) . This has been 
done by subtracting an average background, that is the average sampled during 
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a given energy/mass scan cycle, from the average ion count rates in each energy 
and angle bin. As a consequence, normally non-negative moments such as number 
density may end up negative, when the count rate of a given ion is very close to 
background levels, and the count rate in the background channel happens to be on 
the high side due to normal statistical or temporal fluctuations. This is to be 
expected, and negative values ought to be included in any statistical averaging 
of number densities from these files, in order not to bias the result. 

All MD moments have a standard deviation assigned to them. This one accounts 
for purely statistical uncertainties, those associated with Poisson counting 
statistics. In all cases but one, the tabulated value is an integer number be- 
tween 0 and 999, which represents the ratio in percent (%) between the standard 
deviation and the absolute value of the moment itself, rounded downward (values 
greater than 999 assumed irrelevant) . The one exception is the bulk flow angle 
(drift direction angle) in the GSE X-Y plane (spin plane), where the standard 
deviation is expressed as an angle between 0 and 360 deg, rounded to the nearest 
integer . 

The variance (square of standard deviation) of a given moment, or a given com- 
bination of moments (as in bulk velocity and mean energy) , has been calculated 
in a customary fashion by taking the partial derivative with respect to the 
count rate in each energy and angle bin included in the moment (in both numera- 
tor and denominator, where applicable) , squaring the derivative, multiplying it 
with the variance of the associated count rate, and adding such terms over all 
bins. If the same count rate is used twice, in order to replace a missing sam- 
ple in an adjacent bin (see below), its contribution to the variance is adjusted 
so as to reflect the reduced number of independent samplings. All standard de- 
viations, except the one assigned to bulk flow angle (see below) , include a con- 
tribution from the variance of the background measurement. The reason for the 
exception is that the background measurement, although subtracted from all other 
count rates, is a scalar (single number) that cancels out from the calculation 
of flow angle. 

For various reasons, the data may be statistically insufficient to allow a given 
moment (or combination of moments) to be calculated. If the number density has 
been calculated to be a negative number (background measurement too high) , for 
example, it makes no sense to even attempt to calculate a mean energy. And if 
only a few energy channels have been sampled (due to noisy data) , it makes no 
sense to calculate any of the moments, since the output format presumes a cer- 
tain degree of consistency. Whenever a moment calculation fails, the corre- 
sponding standard deviation is set equal to -1. 


7.1 File Format 

The MD-file has been written in groups of five lines (records), using 
formatted sequential FORTRAN WRITE statements as follows: 

The first of five lines is a title line: 

WRITE (8,230) IYYDDD, JSTART, JSTOP, RX, RY, RZ,RT,DZ, IMLAT, TLOCL, 
*BXMD, BYMD, BZMD, BTMD, BETA, IDBETA, BCTR, MDTCR, 

*ENEMAX, IRATE, RSEY, RSEZ 

230 FORMAT (16,215, 5F6 . 1, 14, F5 . 1, 4F8 . 1 , 1PE 10 . 2 , 13 , 1PE10 .2, 13, 

*0PF5 . 1, 12, 2F6 . 1) 



The variable names represent the following quantities: 

IYYDDD is the year (two digits) and day of year (three digits) . 

JSTART is the universal time in minutes at the beginning of the averaging 
interval, that is the time of the first good data in that interval. 

This is normally at the start of an energy/mass scan cycle, unless 
some initial data in that cycle are bad. 

JSTOP is the universal time in minutes at the end of the averaging interval, 
that is the time of the last good data in that interval. This is 
normally at the end of an energy/mass scan cycle, unless the last 
data in that cycle are bad. 

RX, RY, RZ, and RT are, respectively, the GSM X, Y, Z, and radial distance 
at the midpoint of the averaging interval, all in units of "earth 
radii" ("Re"). RY and RZ are set to 999. if GSM coordinates not 
available (RX same in GSE ) . 

DZ is the distance in "earth radii" (at midpoint) from the nominal neutral 

sheet in the geotail according to Fairfield and Ness [J.Geophys. Res., 
75, 7032, 1970]. This is only displayed for GSM X < -11 Re, otherwise 
set to 999. If no GSM coordinates available, it is set to 0. 

IMLAT is the geomagnetic latitude in degrees (at midpoint) , rounded to the 
nearest integer. This is set to 0, if no ephemeris tape available. 

TLOCL is the geographic local time in hours and 1/10 hours (at midpoint) . 

This is set to 0.0, if no ephemeris tape available. 

BXMD, BYMD, BZMD, and BTMD are the GSE components and absolute value, respec- 
tively, of the measured magnetic field, each averaged over the whole 
MD interval (for data quality check, see next section) . The unit is 
"nanotesla" ("nT"), or equivalently, "gamma". 

BETA and IDBETA are a simplified representation of the ion plasma beta and its 
standard deviation (% of absolute value) . Its definition is explained 
in the next section. 

BCTR is the average background count rate in counts/sec. 

MDTCR is a flag showing which of two detector pulse height triggering levels 

has been used ( = 2 in data from before March 2 of 1980, = 1 afterwards) . 

ENEMAX is the maximum energy sampled, in units of "keV/e", or equal to 16.1, 

whichever is smaller (16.1 is maximum in moments) . Even if it is listed 
as 16.1 (typical) the moments of some ion species may sometimes be lim- 
ited to lower energies, depending on the energy/mass scan mode (see next 
section) . 

IRATE is a flag showing which of two data accumulation (and telemetry trans- 
mission) rates has been used. The low or normal rate (80% of the time) 
is shown by IRATE = 1, the high rate by IRATE * 4. Note: low rate means 
four (4) samplings/sec, high rate means sixteen (16) samplings/sec. 

RSEY and RSEZ are the GSE Y and Z (at midpoint of averaging interval) in units 
of "earth radii". 
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The next four lines list, respectively, the moments for H+ (K= 1) , He++ (K= 2) , 
He+ (K= 3), and 0+ (K= 4) (that is actually for M/Q= 1, 2, 4, and 16) : 


DO 270 K-l, 4 

WRITE (8, 240) DNSS, IDO, DNS 8, ID1,EMN8, ID2, 

*VDRFT, ID3, IDRFT, ID4, 

*BGD, ID6 , ACS EC, XCSEC, 

*IPAMIN, IPAMAX, DENS 8 , ID7,EPER8, ID8,EPAR8, ID9 

240 FORMAT (4 (1PE10 .2, 13) , 15, I 4 , 1PE10 . 2 , 13 , 0PF7 . 2 , 0PF6 . 2 , 

*14,14, 1PE10 .2, 13,2 (1PE9.2, 13) ) 

270 CONTINUE 

The variable names represent the following quantities: 

DNSS and IDO are the number density in "/cm3" and standard deviation (% of ab- 
solute value) of ions with energies between 10 eV/e and about 100 eV/e, 
that is of ions in the lowest energy channel (with RPA voltage fixed at 
10 V) . 

DNS 8 and ID1 are the number density in " /cm3” and standard deviation {% of ab- 
solute value) of ions with energies between about 100 eV/e and 16 keV/e 
(normally) ^ 

EMN8 and ID2 are the mean energy in "keV" and standard deviation (%) of ions in 
the nominal energy range from 100 eV/e to 16 keV/e. Note: this is total 
energy, including bulk motion, and it is in units of "keV" , not M keV/e". 

VDRFT and ID3 are the common (among those ions) drift speed (bulk flow speed) 

in "km/ sec" and standard deviation (%) of ions in that same energy range 
(100 eV/e - 16 keV/e) . This drift speed is in the s/c spin plane, that 
is approximately in the GSE X-Y plane. 

IDRFT and ID4 are the longitudinal direction angle of that same drift and its 

standard deviation, both in "degrees". IDRFT = 0 is sunward and = 90 is 
duskward (-180 deg to +180 deg) . ID4 is between 0 and 360 deg. 

These moments, from DNSS through IDRFT, assume that the velocity distribution 

has rotational symmetry around the drift (flow) vector (in the GSE X-Y plane) . 

BGD and ID6 are the equivalent isotropic density in ”/cm3", over the 100 eV/e 

to 16 keV/e range, and standard deviation (%) corresponding to the meas- 
ured average background count rate. That is, BGD is equal to the total 
background correction of DNS8 (a number already subtracted from DNS8) . 
The ID6, when expressed in absolute terms, is part of ID1 . 

ACSEC is the number of seconds of data included in the moments for that ion 

species. This number is normally a small fraction of the length of the 
averaging interval (energy/mass scan cycle) . 

XCSEC is the number of seconds during the averaging interval when the MD detec- 
tor was in low sensitivity mode and, at the same time, was sampling that 
particular ion species. This number is almost always zero in the magne- 



tosphere, but usually nonzero for H+ in the magnetosheath and solar wind. 
In the latter cases it is often nonzero for He++ ions as well. If it is 
not zero, it implies that at least some of the count rates of that ion 
species have been multiplied by a crude scaling factor (= 65 for MDTCR 
= 2) to compensate for the reduced sensitivity. 

IPAMIN and IPAMAX are the minimum and maximum pitch angles in "degrees" sampled 
for those ions, rounded to the nearest integer. Note: if the magnetic 
field is properly measured and corrected for magnetometer offsets, then 
the sum of these angles, IPAMIN + IPAMAX, should range between about 
170 deg and 190 deg. 

DENS8 and ID7 are the number density in "/cm3" and standard deviation (% of ab- 
solute value) of ions in the 100 eV/e - 16 keV/e range, assuming that the 
velocity distribution has rotational symmetry around the local magnetic 
field vector (see next section for further explanations) . 

EPER8 and ID8 are the mean perpendicular energy (perpendicular to magnetic 

field; two degrees of freedom) in "keV" and standard deviation (%) of 
ions in the 100 eV/e - 16 keV/e range, assuming that same kind of 
symmetry . 

EPAR8 and ID9 are the mean parallel energy (parallel to magnetic field; one 
degree of freedom) in "keV" and standard deviation (%) of ions in the 
100 eV/e - 16 keV/e range, assuming that same kind of symmetry. 


7.2 Further Explanations of Certain Variables. 

The magnetic field values provided here may to some extent duplicate what is 
now available at the NSSDC (Prof. C. T. Russell; NSSDC identification number 
77-102A-0 4 ) and on the World Wide Web (http://www-ssc.igpp.ucla.edu/forms/isee/) 
but are not intended to supplant those. The values derived here may not have 
been adequately screened against telemetry noise or adjusted for instrument 
anomalies, so it is recommended that the PI provided data be consulted whenever 
there is doubt about accuracy. There are two types of errors known to occur 
here : 

a. When the magnetometer is commanded to low gain on the inbound leg of 
the s/c orbit, usually between R * 8 Re and = 5 Re, the corresponding 
flag may show up slightly too late on the raw data tape, resulting in 
a brief underrepresentation of the field by a factor of 32 . This may 
last for about a minute or less (part of one major frame) and affects 
a single MD interval. This error ought to be fairly obvious in the ED 
data, and therefore traceable in the MD data as well (by comparing the 
times) . The corresponding mismatch at high gain command on the out- 
bound leg, that is a sudden 32-fold increase of the field, is usually 
discovered as "unreal" and corrected. 

b. Part of the computing procedure to decipher the raw telemetry data is to 
infer, indirectly, the current magnetometer antenna flip status from the 
uncorrected field values themselves (pertains to the s/c Y and Z coordi- 
nates) . This part fails occasionally, which in turn makes the remaining 
procedure unable to correct the field readings for magnetometer offset, 
resulting in anomalous modulations of the GSE components. Early in the 
mission, when the offset is fairly small, these modulations cause the 
sum BXMD * *2 +BYMD* *2+BZMD * *2 to be much smaller than the square of the 



time-averaged absolute field strength BTMD* *2 . Later in the mission, 
when magnetometer drift has created a larger offset, this sum may settle 
at a rather large artificial value, of order 100 nT . In either case, it 
is recommended that the two measures of field strength be compared rou- 
tinely, and that other magnetic field data be consulted if those two 
differ by more than a few percent. This problem usually only occurs in 
regions of weak magnetic field, mostly in the solar wind, where it is 
easily recognized by the aforementioned sum being nearly constant, with 
no evidence of diamagnetic influence by the ion pressure (see below) . 

An additional (or alternative) check can be made by summing IPAMIN and 
IPAMAX. If that sum is several degrees smaller than 170 deg or several 
degrees larger than 190 deg, then the magnetic field should be in doubt. 

The velocity moments have been calculated two ways, A and B, using either of two 
energy-angle matrices. Both matrices consist of 32 energy channels (covering 
the entire range from 10 eV/e to about 18 keV/e) , but one has 12 spin angle sec- 
tors (30 deg each) and no pitch angles, the other has 9 pitch-angle sectors (20 
deg each) and no spin angles. In method A, using 12 spin angle sectors, the 
calculations consist of the following steps: 

A. I A phase space density and a corresponding standard deviation are 
assigned to each matrix point that has been sampled, using the 
local average count rate, minus an average background count rate, 
and the number of samplings. No interpolations are made at this 
stage . 

A. II Within each energy channel that has been sampled (usually all have), 

except the bottom and top ones, the phase space densities are weighted 
by cosines or sines of the spin angle (center of sector) and summed 
over angle to form two orthogonal projections. At this stage the 
phase space densities are interpolated in angle, if some angular bins 
have no samples. These projections are in turn weighted by energy and 
by an energy bandwidth (see below) and summed over energy channels, 
forming two orthogonal components, approximately the GSE X and Y com- 
ponents, of a vector that is proportional to number flow density. 

Both components are then divided by a total (scalar) sum of phase 
space densities weighted by the energy bandwidth and by the square 
root of the energy, that is by a sum proportional to number density, 
to form (approximate) X and Y drift velocity components, and a drift 
angle IDRFT . Similar summations are carried out with the variances 
(square of standard deviations), using the corresponding partial de- 
rivatives (squared) as weights, to derive the standard deviation of 
the drift angle (ID4) . 

A. Ill A spherical coordinate system is envisioned with its polar axis along 
the drift velocity vector, that is in the GSE X-Y plane. It is now 
assumed that the phase space density has azimuthal symmetry in this 
coordinate system (rotational symmetry around drift velocity vector) . 

The number density, flow density (but not angle) , and energy density 
are recalculated by summing over solid angle and energy in this system. 
The solid angle weighting factors in this case are zones on a unit 
sphere, each defined at the intersection with the GSE X-Y plane by the 
boundaries of a 30 deg spin angle sector, or by one boundary and the 
drift velocity vector. These zones are partially overlapping, and 
summing over 12 spin angle sectors (typically = 14 zones) means cover- 
ing the unit sphere twice, so a factor 1/2 is applied to each sum. 

The drift speed (VDRFT) and mean energy ( EMN8 ) are obtained by dividing 



flow density and energy density, respectively, by the number density 
(DNS 8 ) . No drift energy is subtracted from this mean energy. The cor- 
responding standard deviations (ID3, ID2, and ID1) are derived from 
similar sums of variances, using partial derivatives as weights. 

A. IV The partial number density of ions in the bottom energy channel (DNS5) 

and the corresponding standard deviation (IDO) are also calculated in 
this coordinate system, although in this case the energy "sum” has a 
single term covering the channel width (10 eV/e to about 100 eV/e) . 

If the drift speed in step A. II corresponds to less than 1 eV (less than 14 
km/sec for H+, less than 4 km/s for 0+, etc.), then the coordinate system is 
instead aligned with the X-Y projection of the magnetic field (calculated from 
BXMD and BYMD) , provided the field elevation angle (including BZMD) is between 
-45 and +45 deg. If the latter is not the case, and if the drift speed is be- 
low minimum, the coordinate system is aligned with the s/c spin axis, and the 
phase space densities are treated as isotropic. Only density and mean energy 
are recalculated in these two cases . 

In method B, using 9 pitch-angle ranges (0.0 - 19.9, 20.0 - 39.9, etc.), 
the steps are as follows: 

B. I A phase space density and a corresponding standard deviation are 

assigned to each matrix point that has been sampled, using the 
local average count rate, minus an average background count rate, 
and the number of samplings. 

B.II If only part of the pitch-angle range has been sampled, which is 
often the case (indicated by IPAMIN and IPAMAX) , then the empty 
angular bins near 0 deg and 180 deg are assigned the same phase 
space density as the nearest sampled bin (closer to 90 deg) at 
the same energy. 

B.III A spherical coordinate system is envisioned with its polar axis 

along the magnetic field vector (arbitrary direction) . It is now 
assumed that the phase space density has azimuthal symmetry in this 
coordinate system (rotational symmetry around magnetic field vector) , 
which is to say that the phase space densities in the 9 pitch-angle 
bins represent the entire unit sphere. The number density and the 
parallel (axial) and perpendicular energy densities are calculated 
by summing over solid angle and energy in this system. When summing 
over solid angle, the parallel and perpendicular energies are repre- 
sented by, respectively, the cosine square and the sine square of 
the pitch angle at the center of each bin. The solid angle weighting 
factors in this case are 9 contiguous zones on a unit sphere, each 
defined by the boundaries of a 20 deg pitch-angle sector. The two 
mean energies (EPAR8 and EPER8 ) are obtained by dividing the respec- 
tive energy densities by the number density (DENS8 ) . The correspond- 
ing standard deviations (ID9, ID8, and ID7) are derived from similar 
sums of variances, using partial derivatives as weights. Only the 
standard energy range (0.1 - 16 keV/e) is included here. 

EPER8 and EPAR8 should have a ratio of 2:1 if the velocity distribution is iso- 
tropic (or if the count rates have been extrapolated from a single bin at 90 deg 
pitch angle), since EPER8 represents two degrees of freedom and EPAR8 only one. 
Due to rounding-off errors in the summation over angle, however, the isotropic 
ratio is not exactly 2, but about 1.969. More specifically, if statistical and 
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instrumental errors are neglected, EPER8 is very nearly exact (about 0.02% too 
small), but EPAR8 is about 1.55% too large. Naturally, neither energy will be 
very meaningful if the magnetic field is in error (see above) . 

The "ion plasma beta" listed in the title line, BETA, is actually calculated 
with method A. The ion pressure used here is equal to 2/3 (two degrees of free- 
dom) of the sum of "thermal" energy densities of the four ion species, that is 
total energy densities minus the energy densities associated with the respec- 
tive drift speed, as calculated assuming symmetry around the respective drift 
velocity vector. This pressure has been divided by a magnetic pressure based 
on BTMD . An alternative and more formally correct value may be calculated using 
the "gyrotropic” density DENS8 and perpendicular energy EPER8, although those 
quantities do not account for bulk flow (drift) . In case the magnetic field 
is suspect, and other field measurements are available, the beta can only be 
recalculated with method A (using DNS8, EMN8, and VDRFT for each ion) . 

The summation over energy treats each energy channel, except the first and 
second ones, as a point measurement at the center energy, and takes the energy 
bandwidth to be the distance between adjacent sampled channels in the trape- 
zoidal fashion, ignoring intermediate channels with no samples. At the second 
channel an extra term is added to extend the standard energy integral from the 
center of the channel downward to about 0.1 keV, assuming the flux to be a con- 
stant. This addition brings the mathematical energy range into better agreement 
with the instrumental range of acceptance. At the first energy channel, which 
is treated separately, the energy summation has only one term that includes the 
channel width (10 eV/e - 100 eV/e) as a factor (only DNSS and IDO calculated) . 

Since the summations are done after the completion of an energy/mass scan cycle, 
the energy channels have normally been sampled in a contiguous fashion, but the 
summation procedures are set up to accept gaps of as many as four channels (due 
to noisy data), before declaring the data insufficient. 

Depending on the scan mode, the lowest and highest energy channels sampled may 
vary, and are sometimes different for different ions. The summation procedure 
uses the actual lower and upper channels, so the moments of different ions may 
on occasion cover different energy ranges. There are no flags to separate ions 
in that regard in the MD-file (insufficient space), but the actual energy cov- 
erage is shown for each ion in the MC-file (see below) . See also the explana- 
tion of ENEMAX above. 

The variable ACSEC is a sum of elementary time segments associated with each 
commanded setting of the power supplies controlling energy and mass channels. 
These elementary times are 1/4 sec during low bit rate operation (about 80% of 
the time) and 1/16 sec during high bit rate operation. These times include the 
resetting of the power supplies, however, and are slightly longer than the times 
associated with particle counting. For simplicity, the particle counting is 
interrupted for about 12% of the elementary time segments in both low and high 
bit rate operation to allow for resetting. The ACSEC therefore exceeds the ac- 
tual particle counting time by about 14% in both cases. 


8. THE MC-FILE 

This file is essentially an abbreviated listing of the count rates that have 
been converted to velocity moments in the MD-file. All relevant energy channels 
are represented, but the spin angle bins have been replaced, in each separate 
energy channel, by the peak count, the time of the peak (in seconds), the spin 



angle of the peak (angle of motion of ions) , a measure of the angular width of 
the ion flux, and the spin-averaged count rate with standard deviation. No 
pitch-angle information is listed, but the approximate pitch angle at the peak 
count can be derived from the spin angle and the magnetic field orientation. 

Count rates are listed for five ion species, that is for the same four species 
included in the MD-file: H+, He++, He+, and 0+ (actually for M/Q= 1, 2, 4, and 
16), and, in addition, for doubly charged oxygen: 0++ (M/Q= 8) . The listing for 
0++ is usually incomplete (contains fill) because it has only been measured over 
part of the energy range, or it has not been measured at all, depending on the 
energy/mass scan mode. It is common to limit the measurement of the 0++ to the 
lowest 24 energy channels (to energies below about 10 keV/e) , to avoid contami- 
nation by the more abundant 0+, whose M/Q response function partially overlaps 
with that of the 0++ at the highest energies. 

Usually the MC-file displays count rates from all 32 energy channels, but if 
the energy/mass scan cycle is limited to a lower portion of the energy range 
(as may be the case when the instrument has been commanded for solar wind or 
magnetosheath observations), then the MC-file is also limited to lower ener- 
gies. In such cases, it is not uncommon to have a different energy coverage 
for different ions, but the MC-file lists the same number of channels (same 
number of lines) for all five ion species, displaying dummy fill where data 
are missing. The fill consists of -1 or -1.00. 

Note: The MC-file is the only documentation of the actual energy ranges used 

in the MD-file, in case different ions have been measured over different energy 
ranges. The MD moment calculations will extend across missing intermediate en- 
ergies, up to four contiguous missing channels (see above) , but they will not 
extrapolate below or above the end channels, and they will not extend above 16 
keV/e (or below 100 eV/e in most cases) . The actual ranges can be deduced from 
the MC-file by the lowest and highest channels (lines) that have data (rather 
than fill) . Usually the only places where the energy sweep is shortened and 
made different for different ions are in the solar wind and in the magneto- 
sheath, as predicted from average magnetosphere configurations. 

In order to simplify the file format, all averaged counts are treated as repre- 
senting the whole energy/mass scan cycle, even though different ions were in 
fact sampled at different phases of the cycle (the same ion often more than 
once) . However, the time labels on the peak counts do refer to the respective 
times of sampling, rounded downward to the whole second. 

All MC counts and count rates are raw data, except when the instrument has been 
operating in low sensitivity mode. In that case the peak count and some counts 
included in the spin-averaged rate have been multiplied by a crude scaling fac- 
tor (= 65) to compensate for the reduced sensitivity. Such counts are flagged 
accordingly (see below) . None of the counts or count rates has been corrected 
for background counts due to penetrating radiation (mostly MeV electrons and as- 
sociated bremsstrahlung) . However, the average background count rate measured 
during each energy/mass scan cycle is listed as average counts per sample, along 
with the sampling rate, which is expressed in number of samples per second here. 

The MC counts and count rates may be converted to differential fluxes or phase 
space densities using the instrument parameters listed in a separate section 
below . 


8 . 1 File Format 



The MC-file has been written in groups of up to 34 lines (up to 33 records) , 
using formatted sequential FORTRAN WRITE statements as follows: 

The first line in each group has a single character indicating the format of 
the counts (1 or 2), the second line is a title line that also indicates how 
many lines (energy channels) are to follow in the current group: 

WRITE (11, 243) IFORMT, IYYDDD , JSTART, JSTOP , RX, RY, RZ, RT, RSEY, RSEZ, 

*DZ, IMLAT, TLOCL, BXMD, BYMD, BZMD, BTMD, BETA, IDBETA, BGND, IDBGND, 

*MAXES0 , NSPLPS 

243 FORMAT (12/316, 7F6 . 1 , 1 4 , F5 . 1 , 4F8 . 1 , 1PE10 . 2 , 13, 1PE9.2, 13, 213) 


The variable names represent the following quantities: 

IFORMT is the format flag for the counts, showing either of two formats. 

The standard format (IFORMT - 1) lists the averaged count rates in 
floating point with two decimal places, the other one (IFORMT = 2) 
rounds them to the nearest integer in order to allow more positions 
for the peak counts (always integers) on the same line. 

IYYDDD is the year (two digits) and day of year (three digits) . 

JSTART is the universal time in seconds (not minutes) at the beginning of 
the averaging interval, that is the time of the first good data in 
that interval. This is normally at the start of an energy/mass scan 
cycle, unless some initial data in that cycle are bad. 

JSTOP is the universal time in seconds (not minutes) at the end of the 

averaging interval, that is the time of the last good data in that 
interval. This is normally at the end of an energy/mass scan cycle, 
unless the last data in that cycle are bad. 

RX, RY, RZ, and RT are, respectively, the GSM X, Y, Z, and radial distance 
at the midpoint of the averaging interval, all in units of "earth 
radii” ("Re"). RY and RZ are set to 999. if GSM coordinates not 
available (RX same in GSE) . 

RSEY and RSEZ are the GSE Y and Z (at midpoint of averaging interval) in units 
of "earth radii". 

DZ is the distance in "earth radii" (at midpoint) from the nominal neutral 

sheet in the geotail according to Fairfield and Ness [J.Geophys. Res., 
75, 7032, 1970] . This is only displayed for GSM X < -11 Re, otherwise 
set to 999. If no GSM coordinates available, it is set to 0. 

IMLAT is the geomagnetic latitude in degrees (at midpoint), rounded to the 
nearest integer. This is set to 0, if no ephemeris tape available. 

TLOCL is the geographic local time in hours and 1/10 hours (at midpoint) . 

This is set to 0.0, if no ephemeris tape available. 

BXMD, BYMD, BZMD, and BTMD are the GSE components and absolute value, respec- 
tively, of the measured magnetic field, each averaged over the whole 
interval (for data quality check, see description of MD-file) . The unit 



is "nanotesla” ("nT”), or equiva lent ly, "gamma”. 


BETA and IDBETA are a simplified representation of the ion plasma beta and its 
standard deviation (% of absolute value) . Its definition is explained 
in the description of the MD-file. 

BGND and IDBGND are the average background count rate in counts per sample 
(not per second) and standard deviation (%) . 

MAXESO is the number of lines to follow in the current group (usually = 32, 
sometimes less) . Each line represents one energy channel, always 
starting with the lowest channel (10 eV/e - 100 eV/e) , whether or not 
that channel was sampled, and running through the channels contiguously. 
MAXESO is the number of the highest energy channel sampled at any phase 
of the energy/mass scan cycle. 

NSPLPS is the number of samples per second taken during the current energy/mass 
scan cycle, either 4 (low bit rate) or 16 (high bit rate) . 


The next N - MAXESO lines list counting data in N contiguous energy channels, 
starting at the bottom channel (1= 1), for H+ (J= 1), He++ (J= 2), He+ (J=* 3), 

0+ (J= 4), and C++ (J= 5) (M/Q= 1, 2, 4, 16, and 8): 

DO 248 1=1, MAXESO 

IF { IFORMT . Eg . 1 ) THEN 
WRITE (11, 246) 

* (ACTR(J) , IDCTR(J) , IPCTR(J) , IUTPC(J) , IPDEG(J) , IWDEG(J) , J«l, 5) 

2 46 FORMAT (X, F7 . 2, 13, 14, 15, 14, 14, 4 (F7 . 2, 13, 13, 15, 14, 14) ) 

ELSE 

WRITE (11, 247) 

* (IACTR(J) , IDCTR(J) , IPCTR(J) , IUTPC (J) , IPDEG(J) , IWDEG(J) , J=l, 5) 

2 47 FORMAT (X,l5,l3, 16, 15, 14, 14, 4(15, 13,15,15, 14, 14)) 

END IF 

248 CONTINUE 

The variable names represent the following quantities: 

ACTR is the average (spin-averaged) count rate in floating point represen- 
tation (IFORMT = 1), that is the average number of counts per sample 
(4 or 16 samples per second) . 

IACTR is an integer approximation of ACTR, that is ACTR rounded to the nearest 
integer. This format (IFORMT =2) is used when some ion (usually H+) at 
some energy (usually around one keV) has a value of ACTR that would cause 
the standard format to overflow (IACTR is capped in the rare event that 
it would also overflow) . This same integer representation is used for 
all ions at all energies, regardless of where the overflow would have 
ocurred. This condition is usually associated with a change of instru- 
ment operation to low sensitivity mode, so the largest IACTR values 
usually contain counts that have been "corrected”, that is multiplied 
by a scaling factor (= 65), and are flagged accordingly (by IDCTR) . 


zt 



IDCTR 


i 3 normally the Poisson standard deviation of ACTR (or IACTR) f expressed 
in percent and rounded downward to an integer value. That is, IDCTR 
represents the ratio (in %) between the square root of the sum of counts 
and the sum of counts (sum = ACTR * No of samples) . In case the instru- 
ment has been in low sensitivity mode at any time while a particular ACTR 
(or IACTR) was being averaged, the IDCTR is instead used as a flag and 
set equal to -9. In that case ACTR (or IACTR) contains "corrected" 
counts and should be treated with some care. 

IPCTR is the maximum number of counts during any single sample (about 1/4 :th 

or 1/16 :th of a second) of a given ion at a given energy. If IDCTR= -9, 
then IPCTR probably contains a "correction factor" (= 65) . 

IUTPC is the (approximate) time of maximum count rate (IPCTR), expressed as 
integer number of seconds after the beginning of the current averaging 
interval. That is, the universal time in seconds is = JSTART + IUTPC. 

IPDEG is the angle of motion in the s/c spin plane of those ions that cause 
the maximum count rate. The angle refers to the midpoint of the samp- 
ling interval (and the midpoint of the instrument field of view) , and 
is rounded to the nearest integer. During each sampling interval the 
instrument sweeps (spins) through about 30 deg of angle in low bit rate 
(4 samplings/sec), and about 7.5 deg in high bit rate (16 samplings/sec) . 

IWDEG is a measure of the angular width (nearest integer) of the ion flux dis- 
tribution around the maximum. It is a sum of two angles, one being the 
closest spin angle from IPDEG with a count rate less than 1/3 of IPCTR, 
the other being the most distant spin angle from IPDEG with a count rate 
at least 1/3 of IPCTR. That is, IWDEG is a crude measure of the full 
width at 1/3 of maximum. It is measured after IUTPC if more than 1/2 
s/c spin remains, otherwise set to -1. If IPCTR is more than 3 times 
greater than any subsequent counts, IWDEG is set to 30 deg in low bit 
rate and 8 deg in high bit rate (or -1 if too close to end of interval) . 


8.2 Relevant Instrument Parameters. 

To convert counts per sample, or CTS, to counts per second, or CTRATE, use 
CTRATE = CTS *RATE /0. 21865 

where RATE * 1.0 in low bit rate, and RATE = 4.0 in high bit rate. 

However, the count rates in the MC-file are raw counts and may need to be ad- 
justed for detector degradation (including peak counts and background) , depend- 
ing on which pulse height triggering level was being used. Count rates from 
before March 2 of 1980 have been obtained at the higher of two levels in the MD 
particle detector, at the MDTCR = 2 level. The lower level, MDTCR = 1, used 
afterwards, provides a more nearly one to one detection level, that is one count 
for every ion entering the detector, but it also admits more false counts due 
to penetrating radiation than does the higher level. The reason for using the 
MDTCR = 2 level is to minimize the background. The lower sensitivity to ions 
is not a problem in itself as long as it is well known, since the count rate 
can be adjusted accordingly, but it has slowly declined over time, making it 
necessary to do periodic in-flight calibrations. These have consisted of inter- 
comparing the count rates at the MDTCR = 1 and « 2 levels, which are both part 



of the instrument output, during times of extremely low background. Fortu- 
nately, the MDTCR = 1 level has shown no degradation, and can be used as a 
standard reference. The following table shows the results of these intercom- 
parisons in the lowest and highest energy channels for the four principal ions 
(date refers to beginning of month) . The intercomparisons for 0++ (ION= 5) have 
been less extensive but suggest that the 0+ (ION- 4) ratios can be used for the 
0++ as well. 

c Ratios of MDTCR= 2 count rate to MDTCR 35 1 count rate in flight; 
c 

c in lowest energy channel: 

DATA CMDO 1/ 


c launch 

- Jan 78, 

Jul 78, 

Jan 79, 

Jul 79, 

Jan 80, 

Jul 80 

c ION= 1 
* 

.90, 

.81, 

.60, 

.48, 

.36, 

.24, 

c I0N= 2 
★ 

.95, 

.90, 

.80, 

.70, 

. 60, 

.50, 

C ION= 3 
* 

.95, 

. 85, 

.70, 

.60, 

.50, 

.40, 

C ION= 4 
* 

.95, 

.89, 

.75, 

.65, 

.55, 

.45/ 

c 

c in highest energy 
DATA CMD32 / 
c launch - Jan 78, 

channel : 
Jul 78, 

Jan 79, 

Jul 79, 

Jan 80, 

Jul 80 

C ION= 1 
* 

.94, 

.93, 

.82, 

.72, 

.61, 

.50, 

c ION= 2 
* 

• 97, 

.95, 

.90, 

.85, 

.80, 

.75, 

c ION= 3 
* 

• 97, 

.95, 

.90, 

.85, 

.80, 

.75, 

c ION= 4 
* 

• 97, 

.95, 

.90, 

.85, 

.80, 

.75/ 


Intercomparisons have been made at many energies, and it appears that the ratios 
vary about linearly with energy channel number (1 through 32) . That relation 
has been used when converting count rates to velocity moments in the MD-file, 
that is, the ratios have been interpolated linearly in energy channel number. 
They have also been interpolated linearly in time between the dates above. For 
dates beyond February of 1980 this scheme has been abandoned, and MDTCR = 1 
counts have been used exclusively. Given count rates from before March 2, 1980, 
these are thus to be adjusted by 

CTRATE - CTRATE /CMD(IE, date, ION) 

where IE is the energy channel No, and CMD ( . . . ) is obtained by linear interpo- 
lation in IE and date between the CMDO 1 and CMD32 in the table. For ION= 5 
(0++ ) use the same numbers as for ION= 4 (0+) . 

There is actually one more adjustment that can be made to improve the one to one 
relationship between counts and ions in the MC-file. This one does not depend 
on the MDTCR level, but it depends somewhat on the energy /mass scan mode. To 
have exactly one count for every ion entering the MD detector, on average, re- 
quires that the power supplies controlling the M/Q separation be tuned exactly 
to the peak response for a given ion at every energy. This is impractical, but 



in all scan modes, except the so called "load mode", the single mass channel 
chosen to represent a given ion will have let that ion through at about 90% of 
peak response, or higher. In the "load mode", however, every mass channel has 
been sampled, and all counts at a response of 40% or higher have been used (for 
the five ions listed) . The average response for those counts ranges between 70% 
and 90% of the peak, varying somewhat randomly from one energy channel to the 
next, but is mostly 80% to 85%. 

Therefore, if the counts have been obtained in "load mode" it may be worth ad- 
justing the time averaged counts (ACTR and IACTR) by another factor of about 1.2 
(average over energy) . The peak count rates (IPCTR) cannot be adjusted for this 
effect at this stage, since the mass channel No. is not listed. The "load mode" 
can be recognized from the MS-file (see below) . By contrast, all counts used 
when calculating moments for the MD-file were adjusted for off-peak response on 
a sample by sample basis (after first summing over the mass peak in load mode) . 

Once adjusted, the count rates in the MC-file can be converted to differential 
flux "FLUX" and phase space density "F" with the following subroutine. 


SUBROUTINE AFLUX ( ION, CTRATE, IE, FLUX, F ) 
c**** input: ION- 1 (H+) , 2 (He++) , 3 (He + ) , 4 (0+) , or 5 (0++) 
c**** CTRATE— counts per second (floating point) , and 

c**** IE— 1, 2, 3, 32 (energy channel) 

c * * * * 

c**** Output: "FLUX" in units of "/cm2/sec/ster/keV" 
c**** "F'^in units of "sec3/km6". 

Q**** 

DIMENSION AM (5) , Q ( 5) ! ion mass and charge units 

DATA AM/ 1., 4. ,4. ,16. ,16./, Q/l.,2.,1.,1.,2./ 
c**** instrument energies (center of channels) : 

DIMENSION ENERGY (32) ! "keV/e" 

DATA ENERGY/ 


if 

. 040, 

.212, 

.410, 

. 628, 

.851, 

1.095, 

1.353, 

1.633, 

★ 

1.929, 

2.244, 

2.580, 

2.934, 

3.317, 

3.718, 

4.146, 

4.599, 

if 

5.080, 

5.592, 

6.132, 

6.713, 

7.333, 

7.998, 

8.701, 

9 . 446, 

if 

10.235, 

11.076, 

11.969, 

12 . 917, 

13.927, 

14 . 999, 

16.144, 

17.364/ 


c**** instrument geometric factors, including delta-E: 

DIMENSION G (32, 5) ! "1.0E-4 cm2 keV" 

DATA G/ 
c**** ION= 1: 

* 3. 60, 6. 00, 6. 00, 6. 00, 6. 01, 6. 01, 6. 01, 6. 01, 6. 01, 6. 01, 6. 02, 6. 02, 6. 02, 

* 6.48,6.94,7.40,7.85,8.31,8.77,9.23,9.69,10.3,10.9,11.5,12.1,12.8, 

* 13.4,14.0,14.6,15.2,15.8,16.4, 
c**** ION= 2: 

* 4. 50, 7. 37, 7. 20, 7. 04, 6. 87, 6. 70, 6. 54, 6. 37, 6. 21, 6. 04, 5. 87, 5. 71, 5. 54, 

* 5.85,6.17,6.48,6.80,7.11,7.43,7.74,8.06,8.39,8.71,9.04,9.37,9.69, 

* 10.0,10.4,10.8,11.2,11.6,12.0, 
c**** ION= 3: 

* 7. 20, 11. 9, 11. 3, 10. 7, 10. 0,9. 41, 8. 79, 8. 16, 7. 54, 6. 92, 6. 30, 5. 67, 5. 05, 

* 5.37,5.69,6.01,6.32,6.64,6.96,7.28,7.60,7.99,8.38,8.77,9.16,9.54, 

* 9.93,10.3,10.7,11.1,11.5,11.9, 
c**** ION= 4: 

* 3.10,5.28,5.23,5.17,5.12,5.07,5.01,4.96,4.91,4.86,4.80,4.75,4.70, 

* 4.85,4.99,5.14,5.29,5.44,5.59,5.73,5.88,6.26,6.63,7.02,7.42,7.82, 

* 8.23,8.61,9.02,9.42,9.82,10.2, 
c**** ION= 5: 


l-f 



* 4. 71,7. 83,7. 55, 7. 28,7. 01, 6. 74, 6. 47, 6. 19, 5. 92, 5. 65, 5. 38,5. 10, 4. 83, 

* 5.05,5.26,5.47,5.69,5.90,6.12,6.33,6.54,6.92,7.30,7.68,8.06,8.52, 

* 8.96,9.34,9.72,10.2,10.6,11.0/ 
c**** normalize: 

GDE= G ( IE, ION) *1 . OE-4 
c**** differential flux: 

FLUX- CTRATE/GDE/Q (ION) ! " /cm2 / s /sr /keV) " 
c**** phase space density: 

F= 0.5449208 * FLUX* AM ( ION) *AM ( ION) / Q ( ION) /ENERGY (IE) ! "s3/km6" 

RETURN 

END 


The energy channel widths (in keV/e) and the angular fields of view in each 
energy channel are inherent in the geometric factors listed in this subroutine 
(G-delta-E) , but it may be of interest to know them separately: 

The external energy bandwidth is defined by the internal energy resolution of 
the instrument, which is about a constant 5% at all energies, but the external 
bandwidth is not a constant fraction of energy, because the ions are pre-accel- 
erated by about 3.0 kV before they enter energy analysis (energy selection) . 

To obtain the absolute energy bandwidth to incoming ions, add 3.0 keV/e to the 
center energies listed in the subroutine, except the lowest energy, and take 5% 
of the sum. The lowest energy channel is different, due to the applied RPA vol- 
tage, and extends approximately between 10 eV/e and 100 eV/e (more precisely to 
about 110 eV/e) . 

* 

As mentioned, the center of the instrument field of view (the one used here) 
points 5 deg below the spin plane, that is about 5 deg below the GSE X-Y plane 
(ion velocity vector pointing 5 deg above) . In this plane the width is about 
10 deg. In the perpendicular (GSE Z) direction it varies with energy (due to 
the pre-acceleration), from about 45 deg at 10 eV/e to 10 deg at 18 keV/e. To 
be more specific, if the center energies of the 32 energy channels are used as 
reference, the corresponding 32 angular widths are as follows (full width at 
25% of max; in deg) : 


o 

o 

CO 

o 

o 

26.0, 

22.5, 

20.0, 

18.0, 

16.5, 

15.5 

15.0, 

14 .0, 

13.5, 

13.0, 

13.0, 

12 .5, 

12.5, 

12.0 

12.0, 

11.5, 

11.5, 

11.0, 

11.0, 

11.0, 

11.0, 

10.5 

10.5, 

10.5, 

10.5, 

10.5, 

10.5, 

10.0, 

10.0, 

10.0 
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This file contains complete mass spectra in four energy ranges, provided the 
instrument energy/mass scan included complete mass scans at some energies (not 
always the case) . The counts in this file are also from the MD detector, but 
they are normally not the same counts as those listed in the MC-file, which went 
into moment calculations for the MD-file. It is only in the mode called "load 
mode" that the same raw counts have been used in all three files with MD data. 
This particular mode (hardwired) consists of one complete 32-step energy scan in 
each successive mass channel over the complete 64 mass-channel range (four such 
scan cycles have been used to define a full cycle in high bit rate here) . This 
scan mode is used infrequently, and can be recognized by the fact that the num- 
ber of samplings in the MS-file (see below) add up to cover the entire averaging 
interval (entire scan cycle) . In all other scan modes the mass spectra are ob- 
tained as separate and small parts of each energy/mass scan cycle, and are only 


is- 



used to provide data for the MS-file here. For simplicity, the same beginning 
and ending times are listed here as are listed in the MD-file, even if the mass 
spectra were obtained for brief intervals well inside those times. 

The four mass spectra are listed in four vertical columns, with each line repre- 
senting one mass channel, beginning with channel No. 1 at the top and ending 
with No. 64 at the bottom. Each column represents one energy range, with the 
energy increasing from left to right: 

Leftmost energy range: channel No. 1 (about 10 eV/e - 100 eV/e) . 

Second to the right: channels No. 2-5 (about 100 eV/e - 1 keV/e) 

Third to the right: channels No. 6-12 (about 1 keV/e - 3 keV/e) 

Fourth to the right: channels No. 13 - 31 (about 3 keV/e - 16 keV/e) 

The counts in each mass channel (on each line) are a sum of counts from all 
energy channels sampled within each energy range. Along with the sum of counts 
is the number of samples. 

The MS-counts are pure raw counts; no "corrective” factor has been applied when 
the instrument has been operating in low sensitivity mode. However, each column 
has a flag above it showing how many of the samples in that column were taken in 
low sensitivity mode (without specifying mass channel) . No background counts 
have been subtracted either. 

Within each mass spectrum (each column) certain mass peaks may be found more 
often than others. The following is a listing of expected mass channel No. 

(1 through 64) at peak count rate for certain important M/Q: 

M/Q= 16.00 about 16 to 17 

= 8.00 about 21 to 22 

= 4.00 about 29 to 30 

= 2.66 about 35 to 36 

= 2.00 about 40 to 41 

= 1.50 about 45 to 46 

= 1.00 about 55 to 56 


9.1 File Format 

The MS-file has been written in groups of 68 lines (65 records), 
using formatted sequential FORTRAN WRITE statements as follows: 

The first line in each group is blank, the second and third make up a title, 
and the fourth lists (above each mass spectrum) the No. of samples taken in 
low sensitivity mode: 

WRITE (12, 255) 

*IYYDDD, JSTART, JSTOP, RX, RY, RZ, RT, RSEY, RSEZ, DZ, I ML AT , TLOCL, 





*BXMD, BYMD, BZMD, BTMD , BETA, IDBETA, BGND, IDBGND, MAXEMS, NSPLPS, 

* (IMSLOW(I) , 1=1, 4) 

255 FORMAT {X/I6, 215, 7F6. 1, 1 4 , F5 . 1/X, 4F8 . 1 , 1PE 1 0 . 2 , 13 , 1PE9.2, 13,213/ 

*X, 4113) 

The variable names represent the following quantities: 

IYYDDD is the year (two digits) and day of year (three digits) . 

JSTART is the universal time in minutes at the beginning of the averaging 
interval, that is the time of the first good data in that interval. 

This is normally at the start of an energy/mass scan cycle, unless 
some initial data in that cycle are bad. 

JSTOP is the universal time in minutes at the end of the averaging interval, 
that is the time of the last good data in that interval. This is nor- 
mally at the end of an energy/mass scan cycle, unless the last data in 
that cycle are bad. 

RX, RY, RZ, and RT are, respectively, the GSM X, Y, Z, and radial distance 
at the midpoint of the averaging interval, all in units of "earth 
radii" ("Re"). RY and RZ are set to 999. if GSM coordinates not 
available (RX same in GSE) . 

RSEY and RSEZ are the GSE Y and Z (at midpoint of averaging interval) in units 
of "earth fadii". 

DZ is the distance in "earth radii" (at midpoint) from the nominal neutral 

sheet in the geotail according to Fairfield and Ness [J.Geophys. Res., 
75, 7032, 1970] . This is only displayed for GSM X < -11 Re, otherwise 
set to 999. If no GSM coordinates available, it is set to 0. 

IMLAT is the geomagnetic latitude in degrees (at midpoint), rounded to the 
nearest integer. This is set to 0, if no ephemeris tape available. 

TLOCL is the geographic local time in hours and 1/10 hours (at midpoint) . 

This is set to 0.0, if no ephemeris tape available. 

BXMD, BYMD, BZMD, and BTMD are the GSE components and absolute value, respec- 
tively, of the measured magnetic field, each averaged over the whole 
interval (for data quality check, see description of MD-file) . The unit 
is "nanotesla" ("nT"), or equivalently, "gamma". 

BETA and IDBETA are a simplified representation of the ion plasma beta and its 
standard deviation (% of absolute value) . Its definition is explained 
in the description of the MD-file. 

BGND and IDBGND are the average background count rate in counts per sample 
(not per second) and standard deviation (%) . 

MAXEMS is the number of the highest energy channel included in the mass spectra 
(31 or lower) . 

NSPLPS is the number of samples per second taken during the current energy/mass 
scan cycle, either 4 (low bit rate) or 16 (high bit rate) . 
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IMSLOW is the total number of samples taken in low sensitivity mode while 
gathering the mass spectrum below it (regardless of mass channel) . 


The next 64 lines, one for each mass channel, list the accumulated number of 
counts (sums) in each of four energy ranges (explained above), along with the 
corresponding number of samples : 

DO 258 J=l, 64 

WRITE (12, 257) ( (ISUM04 (I) , ISPL04 (I) ) , 1*1, 4) 

257 FORMAT (X, 4 (19, 14) ) 

258 CONTINUE 

The variable names represent the following quantities: 

ISUM04 is the accumulated number of counts within energy range I, regardless 
of instrument look angle. 

ISPL04 is the corresponding number of samples (4 per seconds in low bit rate, 

16 per seconds in high; cf . NSPLPS) 

10. THE EX-FILE 

* 

This file contains auxiliary ISEE-1 ephemeris that may be used for labeling 
plots of data from the other files (see comments on existing plots at the end 
of this data user's guide) . It is derived from separate ISEE ephemeris tapes 
which have a more complete set of coordinates, all listed once per minute. Thi 
file has a subset of the ephemeris, and it only lists this set every 30 minutes 
If the ephemeris tape is unavailable (unreadable) , it substitutes the GSE coor- 
dinates from the instrument data telemetry tape for the corresponding ephemeris 
coordinates, choosing these to be as closely as possible on the half hour, and 
fills the remaining positions with dummy numbers. 


10.1 File Format 

The EX-file always has 49 lines (records) which have been written with 
a formatted sequential FORTRAN WRITE as follows: 

MMMM~0 

DO 38 1-1, 49 

36 WRITE (7,37) IYYDDD, MMMM, RX, RY, RZ, RT, DZ, IMLAT, TLOCL, RSEY, RSEZ , RL 

37 FORMAT (16, 15, 5F6 . 1, 14, F5 . 1, 2X, 3F6 . 1) 

38 MMMM=MMMM+ 3 0 

The variable names represent the following quantities: 

IYYDDD is the year (two digits) and day of year (three digits) . 

MMMM is the universal time in minutes. 



RX, RY, RZ, and RT are the GSM X, Y, Z, and radial distance, respectively, 
all in units of "earth radii" ("Re") . RY and RZ are set to 999. if 
GSM coordinates not available (RX same in GSE) . 

DZ is the (calculated) distance in "earth radii" from the nominal neutral 

sheet in the geotail according to Fairfield and Ness [J.Geophys. Res., 

75, 7032, 1970] . This is only displayed for GSM X < -11 Re, otherwise 
set to 999. If no GSM coordinates available, it is set to 0. 

IMLAT is the geomagnetic latitude in degrees, rounded to the nearest integer. 
This is set to 0, if no ephemeris tape available. 

TLOCL is the geographic local time in hours and 1/10 hours. This is set to 0, 
if no ephemeris tape available. 

RSEY and RSEZ are the GSE Y and Z in "earth radii" (may have values even when 
no ephemeris tape available; same with RX and RT) . 

RL is Mcllwain's L value. This is only defined in the inner magnetosphere, 

otherwise set to 0. (also when no ephemeris tape available). 

11. GRAPHIC REPRESENTATION 

A subset of the data in the ED- and MD-files is available in graphical form in 
a series of booklets, one booklet for each month of data. There are two pages 
of graphs for eac]} day, each page covering 12 hours of universal time. Each 
page consists of a vertical stack of 7 panels. These panels show the following 
data items : 

Panel 1 (top) shows the "total" ion density from the ED— file (EDDNS) and the 
partial densities of He++, He+, and 0+ from the MD-file (DNSS) . Note: the 
labels He+ and 0+ may not be appropriate in the magnetosheath and solar wind, 
where heavier ions in high charge states may dominate weak count rates near 
those M/Q values, especially near M/Q= 4. 

Panel 2 shows the He++/H+ density ratio from the MD-file (DNS8 [He++] /DNS8 [H+] ) . 
This panel and the top one are the only ones with MD-data, the ones below all 
have data from the ED-file only. 

Panel 3 shows the ion "thermal" energy from the ED-file (EDEMN) . 

Panel 4 shows the "thermal" energy density from the ED-file multiplied by 2/3 
(EDDNS*EDEMN*2 . /3 . ) along with the magnetic pressure and field strength based 
on the ED-averages . Note: the magnetic pressure and field strength are based, 
respectively, on the sum and the square root of BXED**2 + BYED* *2 + BZED**2 
(not on BTED) in order to show more clearly those (rare) times when the magnetic 
field is flawed. When the problem is the one related to the wrong antenna flip 
status, the magnetic pressure typically stays nearly constant, showing no evi- 
dence of the diamagnetic effect of the ions (total pressure is not preserved) . 
This constant value is small in the early data, while the magnetometer offset is 
still small, but increases over time, due to increasing offset. The brief (one 
data point) error in the magnetic field gain factor on many inbound (and a few 
outbound) passes is also clearly seen in this panel (errors discussed in connec- 
tion with the ED- and MD-files above) . 

Panel 5 shows the ion drift speed (in GSE X-Y plane) from the ED-file (EDVEL) . 



Panel 6 shows the magnetic field elevation angle (relative to GSE X-Y plane) 
from the ED-file ( ITHETB) . Note: when the magnetic field is in error, because 
of the antenna flip status being misinterpreted, this angle typically stays 
close to either plus or minus 90 degrees. 

Panel 7 (bottom) shows the ion drift angle ( IEDANG) and magnetic longitude angle 
( IPHIB) in the GSE X-Y plane from the ED-file. 

Below the bottom panel are the UT and assorted ephemeris . When ephemeris tape 
has not been available (very rare) only the geocentric distance (R) will show. 
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A2. Magnetosphere /Geomagnet ically Disturbed Day 
(High CW-/H+ density ratio) 
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Appendix 4 


VAX/VMS Command Files 

for Magnetic Tape Operations 

(13 pages) 



$!************** LOCKHEED PLASMA COMPOSITION EXPERIMENT ON ISEE-1 ************** 
$ ! ****** COPY_ISEEMS.COM ****** 

$! 

$! Interactive command file for copying data from archival tapes (label: ISEEMS) 
$ ! using the VAX/VMS Backup Utility 

$ ! 

$! by O. W. Lennartsson 

$! Lockheed Martin Missiles & Space 

$ ! Advanced Technology Center 

$ ! Dept. HI - 1 1 , Bldg. 255 

$! 3251 Hanover Street 

$ | Palo Alto, CA 94304 

$! 

$! Telephone (650) 424-3259 

$! 

$! Developed under NASA Contract NAS5-33047 and NASA Grant NAGW-4496 

$! 

$! Operating system used: VAX/VMS version V5 . 5 

$ ! 

$! To execute this file, simply enter 0COPY_ISEEMS, if the file resides in the 
$! same directory used for the data files, otherwise specify directory name in 
$! the command: @directory_nameCOPY_ISEEMS (your default directory should be the 
$! one used for data files) . 

$ i ****************************************************************************** 

$ WRITE SYS$OUTPUT " " 

$DEFAULT_DIR = F$LOGICAL ( "SYS$DISK" ) +F$DIRECTORY ( ) 

$ WRITE SYS$OUTPUT„"Your default directory is: ", " DEFAULT_DIR' 

$WRITE SYS$OUTPUT "Is this where you want the files?" 

$ WRITE SYS$OUTPUT "If not, Ctrl/Y, reset default and reenter command." 

$ WRITE SYS$OUTPUT " " 

$INQUIRE DUM "If directory OK, then press return" 

$SET NOON 

$ WRITE SYS$OUTPUT " " 

$INQUIRE NAME "Which tape drive? (enter name of device) " 

$DRIVE= F$STRING (NAME) 

$ ALLOCATE 'DRIVE' TAPE 

$IF $SEVERITY .EQ. 1 THEN GOTO DRIVE_OK 
$ WRITE SYS$OUTPUT "EXIT" 

$WRITE SYS$OUTPUT "Reenter command if trying another drive." 

$SET ON 
$EXIT 
$DRIVE_OK : 

$SET ON 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYS$OUTPUT "Mount tape on ' ' F$LOGICAL ( "TAPE" ) ' " 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE DUM "When tape mounted and on line (!), then press return" 
$MOUNT/FOREIGN/DENSITY=6250 TAPE 
$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT - 

"If you only want to copy the Data User's Guide, enter Y (yes) here," 

$ INQUIRE GUIDE - 
"otherwise press return again" 

$IF GUIDE THEN SET NOON 

$IF GUIDE THEN GOTO ONLY_GUIDE 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT "Data files are copied by date, starting with the 


/ 



"first date entered, " 

$WRITE SYS$OUTPUT "and ending with the second date (inclusive) 

"-assuming one and the same year." 

$ WRITE SYSSOUTPUT "If data run into second year, reenter a new set ",- 
"of dates for that second year." 

$ WRITE SYS$OUTPUT - 

"It is OK to copy one single day also, that is second date = first date." 
$ IRWND=0 
$NYOLD=0 
$NEW_YEAR : 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P2 "What YEAR? (two digits)" 

$LENGTH_OF_YEAR : 

$LP2=F$LENGTH(P2) 

$IF (LP2 . EQ . 2 ) THEN GOTO CHECK_YEAR 

$INQUIRE P2 "Enter two digits for the year: 77 through 82” 

$GOTO LENGTH_OF_YEAR 
$WRONG_YEAR : 

$ INQUIRE P2 "Enter either of 77, 78, 79, 80, 81, or 82" 

$CHECK_YEAR : 

$NY=F$ INTEGER (F$EXTRACT (0,2,P2) ) 

$IF (NY. LT . 77 . OR. NY . GT . 82 ) THEN GOTO WRONG_YEAR 

$IF (NY . EQ . NYOLD) THEN GOTO YEAR_OK 

$ IF (NY . LT . NYOLD ) THEN GOTO WRONG_ORDER_YEARS 

$MNDAY=999 

$MXDAY=0 

$YEAR_OK : 

$IRWRS=0 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P3 "First DAY to be copied? (one to three digits)" 

$LP3=F$LENGTH (P3) 

$IF (LP3 .GT . 3) THEN GOTO YEAR_OK 
$IF (LP3.EQ.2) THEN P3:= ”0"P3'" 

$IF (LP3.EQ. 1) THEN P3:= "00"P3' n 
$ND=F $ INTEGER ( F$EXTRACT ( 0 , 3,P3) ) 

$IF (NY.EQ. 80) .AND. (ND . LT . 1 . OR . ND . GT . 366 ) .OR. - 

(NY.NE.80) .AND. (ND . LT . 1 . OR . ND . GT . 3 65 ) THEN GOTO YEAR_OK 
$IF (ND.LE.MXDAY) THEN IRWRS=1 
$IF ( IRWRS . EQ . 1 ) THEN GOTO OVERLAP 1 
$NO_OVERLAP 1 : 

$ISTART=ND 
$ YEAR_ST I LL_OK : 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P3 "Last DAY to be copied? (in that same segment)" 

$LP3=F$LENGTH (P3) 

$IF (LP3 . GT . 3 ) THEN GOTO YEAR_STILL_OK 
$IF (LP3 . EQ . 2 ) THEN P3:= ”0"P3'" 

$IF (LP3 . EQ . 1 ) THEN P3:= "00"P3' n 
$ND=F$INTEGER (F$EXTRACT (0, 3, P3) ) 

$IF (NY . EQ . 80 ) .AND. (ND . LT . 1 . OR . ND . GT . 3 6 6 ) .OR. - 

(NY.NE.80) .AND. (ND . LT . 1 . OR . ND . GT . 365 ) THEN GOTO YE AR_S TIL L_0 K 
$IF (ND . LT . I START) THEN GOTO WRONG_ORDER_DAYS 
$IF (ND.GE.MNDAY) .AND. (IRWRS. EQ.l) THEN GOTO OVERLAP2 
$NO_OVERLAP 2 : 

$IRWRS=0 
$ISTOP=ND 
$IDAY=I START 


2 . 


$SET NOON 
$ BACKUP_TO_D I S K : 

$IF ( IDAY . LT . MNDAY) THEN MNDAY=IDAY 
$ IF ( XDAY . GT . MXDAY ) THEN MXDAY=IDAY 

$IF (IDAY . LE . 9) THEN YYDDD= F$STRING (NY) + "00” + F$STRING ( IDAY) 

$IF (IDAY.GT . 9) .AND. (IDAY. LE. 99) THEN - 
YYDDD= F$STRING (NY) + "0" + F$STRING ( IDAY) 

$IF (IDAY.GT . 99) THEN YYDDD= F$STRING (NY) + F$STRING ( IDAY) 

$IF (IDAY.EQ. ISTART) THEN YYDDD1=YYDDD 

$IF ( IRWND . EQ . 0 ) THEN BACKUP/LOG TAPE : ' YYDDD ' . BCK ' DEF AULT_D I R ' 

$IF (IRWND.EQ. 1) THEN BACKUP /REWIND /LOG TAPE YYDDD' . BCK ' DEFAULT_DIR' 

$IRWND=0 

$IDAY=IDAY+1 

$IF ( IDAY. LE . I STOP) THEN GOTO BACKUP_TO_DISK 
$NYOLD=NY 

$ WRITE SYS$OUTPUT " " 

$IF (ISTOP .EQ. ISTART) THEN WRITE SYSSOUTPUT - 
"Last segment copied: a single day ", 'YYDDD' 

$IF (ISTOP .GT . ISTART) THEN WRITE SYS$OUTPUT - 
"Last segment copied: days ", ' YYDDD1' , " through ",' YYDDD' 

$ WRITE SYSSOUTPUT ” ” 

$ INQUIRE MORE - 

"Are there more data files to be copied from this tape? [Y/N]" 

$IF MORE THEN GOTO NEW_YEAR 
$ WRITE SYSSOUTPUT ” " 

$ INQUIRE GUIDE - 

"Do you want to copy the User's Guide? (last file on the tape) [Y/N] 
$ONLY_GUIDE : 

$IF GUIDE THEN - 

BACKUP/LOG TAPE: GUIDE. BCK ' DEFAULT_DIR' 

$IF SSEVERITY . EQ . 1 THEN GOTO REMO VE_T AP E 

$DIR 'DEFAULT_DIR' GUIDE. 

$IF SSEVERITY . EQ . 1 THEN GOTO REMOVE_T AP E 
$IF GUIDE THEN WRITE SYSSOUTPUT " " 

$IF GUIDE THEN WRITE SYSSOUTPUT "Tape already at end, must be rewound 
$IF GUIDE THEN - 

BACKUP /REWIND /LOG TAPE : GUIDE . BCK ' DEFAULT_DIR' 

SGOTO REMO VE_T APE 
$ WRONG_ORDER_D AY S : 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "First and last days in wrong order!" 

SINQUIRE REDO "Want to reenter dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
SGOTO REMO VE_T APE 
$WRONG_ORDER_YEARS : 

S WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "Going backwards in years (rewinding tape)?" 
SINQUIRE SURE "Is this really the year you want? [Y/N]" 

$IF SURE THEN IRWND=1 
$IF SURE THEN NYOLD=0 
$IF SURE THEN MNDAY=999 
$IF SURE THEN MXDAY=0 
$IF SURE THEN GOTO YEAR_OK 

SINQUIRE REDO "Want to correct year? [Y/N]” 

$IF REDO THEN GOTO NEW_YEAR 
SGOTO REMOVE_TAPE 
SOVERLAP1 : 



$ WRITE SYS$OUTPUT " ” 

$ WRITE SYSSOUTPUT "Caution:” 

$IF (ND . LT . MXDAY) THEN WRITE SYSSOUTPUT - 
"First day precedes days already copied!" 

$ IF (ND.EQ. MXDAY) THEN WRITE SYS$OUTPUT - 
"First day has already been copied once!" 

$! IF (ND .EQ. MXDAY) THEN GOTO NEW_YEAR 

$WRITE SYSSOUTPUT "Tape will have to be rewound." 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE SURE "Is that what you want? [Y/N] " 

$IF SURE THEN IRWND=1 

$IF SURE THEN GOTO N0_0VERLAP1 

SINQUIRE REDO "Want to correct dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
$GOTO REMO VE_T APE 
SOVERLAP2 : 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "Caution:" 

SIF(ND.GT.MNDAY) THEN WRITE SYSSOUTPUT - 
"Check your Last day entry!" 

$IF (ND.EQ. MNDAY) THEN WRITE SYSSOUTPUT - 
"Last day has already been copied once!” 

$! IF (ND.EQ. MNDAY) THEN GOTO YE AR_S T I LL_OK 
$ WRITE SYSSOUTPUT " " 

SINQUIRE SURE "Are you satisfied with your dates? [Y/N]" 

$IF SURE THEN GOTO NO_OVERLAP2 

SINQUIRE REDO "Waot to reenter both dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 

SINQUIRE REDO "Want to reenter last day? [Y/N]" 

$IF REDO THEN GOTO YEAR_STILL_OK 
$ REMO VE_T APE : 

SSET ON 

SDISMOUNT 'DRIVE' 

SDEALLOCATE 'DRIVE' 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "EXIT" 

$ WRITE SYSSOUTPUT "That's the end of this run." 

$ WRITE SYSSOUTPUT - 

"Tape is dismounted and ' ' FSLOGICAL ("TAPE" ) ' is deallocated. 
SEXIT 


V 



************** LOCKHEED PLASMA COMPOSITION EXPERIMENT ON ISEE-1 ************** 

****** USER_TAPE.COM ****** 

Interactive command file for copying ISEEMS data files from disk to a new 
magnetic tape using the VAX/VMS Copy Utility 

by 0. W. Lennartsson 
Lockheed Martin Missiles & Space 
Advanced Technology Center 
Dept. HI— 11/ Bldg. 255 
3251 Hanover Street 
Palo Alto, CA 94304 

Telephone (650) 424-3259 

Developed under NASA Contract NAS5-33047 and NASA Grant NAGW-4496 
Operating system used: VAX/VMS version V5 . 5 


To execute this file, simply enter @USER_TAPE, if the file resides in the 
same directory used for the data files, otherwise specify directory name in 
the command: @directory_nameUSER_TAPE (your default directory should be the 
one used for data files) . The new tape will be labeled IONS and will be 
written at a density of 1600 bpi . 

****************************************************************************** 
$ WRITE SYS$OUTPUT " " 

$DEFAULT_DIR = FSLOGICAL ( "SYS$DISK" ) +FSDIRECTORY ( ) 

$ WRITE SYS$OUTPUT "Your default directory is: " , ' ' DEFAULT_DIR' 

$ WRITE SYSSOUTPUT "Is this where the ISEEMS files are?" 

$ WRITE SYSSOUTPUT "If not, Ctrl/Y, reset default and reenter command." 

$ WRITE SYSSOUTPUT " " 

SINQUIRE DUM "If directory OK, then press return" 

$SET NOON 

$ WRITE SYS$OUTPUT " " 

SINQUIRE NAME "Which tape drive? (enter name of device)" 

$DRIVE= FSSTRING (NAME) 

$ ALLOCATE 'DRIVE' TAPE 

$IF SSEVERITY .EQ. 1 THEN GOTO DRIVE_OK 
$ WRITE SYSSOUTPUT "EXIT" 

$WRITE SYS$OUTPUT "Reenter command if trying another drive." 

$SET ON 
$EXIT 
$DRIVE_OK : 

$SET ON 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "Mount tape on ' ' F$LOGICAL ( "TAPE" ) ' " 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT "Make sure tape is on line." 

$ WRITE SYSSOUTPUT ” " 

SINQUIRE INIZLD - 


"Has tape already been initialized at 1600 bpi and labeled IONS? [Y/N] 
$IF INIZLD THEN GOTO INITIALIZED 

$INITIALIZE/OVERRIDE=EXPIRATION/DENS = l 600 / PROT= (G:RW, W:RW) - 
'DRIVE' IONS 
SINITIALIZED: 

$MOUNT/DENSITY=1600 'DRIVE' IONS TAPE 
$ WRITE SYSSOUTPUT " " 
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$ INQUIRE GUIDE - 

"Do you want to copy the Data User's Guide at this point? [Y/N]" 

$ IF .NOT. GUIDE THEN GOTO DATA_NEXT 
$SET NOON 

$COPY/LOG 'DEFAULT_DIR' GUIDE. TAPE 

$IF $SEVERITY .EQ. 1 THEN GOTO DATA_FILES 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE DATA "Do you want to copy data files anyway? [Y/N]" 

$ IF .NOT. DATA THEN GOTO TAPEDIR_LIS 
$SET ON 

$GOTO DATA_NEXT 
$DATA_FILES: 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE DATA "And you also want to copy data files? [Y/N]" 

$IF .NOT. DATA THEN GOTO TAPEDIR_LIS 
$SET ON 
$DATA_NEXT : 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "Data files are copied by date, starting with the " 
"first date entered, " 

$WRITE SYS$OUTPUT "and ending with the second date (inclusive) ",- 
"-assuming one and the same year." 

$ WRITE SYS$OUTPUT "If data run into second year, reenter a new set ", 
"of dates for that second year." 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT - 

"If you are copying one single day, then make last date = first date 
$NYOLD=0 
$NEW_YEAR : 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P2 "What YEAR? (two digits)" 

$LENGTH_OF_YEAR : 

$LP2=F$LENGTH (P2) 

$IF (LP2 .EQ.2) THEN GOTO CHECK_YEAR 

$INQUIRE P2 "Enter two digits for the year: 77 through 82" 

$GOTO LENGTH_OF_YEAR 
$WRONG_YEAR: 

$ INQUIRE P2 "Enter either of 77, 78, 79, 80, 81, or 82" 

$CHECK_YEAR : 

$NY=F$ INTEGER (F$EXTRACT ( 0 , 2 , P2 ) ) 

$IF (NY.LT.77 . OR . NY . GT . 82 ) THEN GOTO WRONG_YEAR 
$IF (NY.EQ.NYOLD) THEN GOTO YEAR_OK 
$YY= F$STRING (NY) 

$DIR/OUTPUT=DRECTRY . TMP ' DEFAULT_DIR' EX' YY' * . DAT 
$IF $SEVERITY .EQ. 1 THEN GOTO YEAR_POSSIBLE 
$DELETE DRECTRY . TMP ; * 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "No files found from year ”,'YY',"!" 

$ INQUIRE REDO "Want to reenter year? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
$GOTO REMO VE_T AP E 
$YEAR_POSSIBLE : 

$DELETE DRECTRY . TMP ; * 

$ IF (NY . LT . NYOLD ) THEN GOTO WRONG_ORDER_YEARS 

$MNDAY=999 

$MXDAY=0 

$YEAR OK: 



$ IRWRS=0 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P3 "First DAY to be copied? (one to three digits)" 
$LP3=F$LENGTH (P3) 

$IF (LP3 .GT.3) THEN GOTO YEAR_OK 
$IF (LP3 . EQ . 2 ) THEN P3:= "0"P3'" 

$IF (LP3 . EQ . 1 ) THEN P3:= "00"P3' n 
$ND=F$INTEGER(F$EXTRACT(0,3,P3) ) 

$IF (NY.EQ. 80) .AND. (ND . LT . 1 . OR . ND . GT . 3 6 6 ) .OR. - 

(NY .NE . 80 ) .AND. (ND . LT . 1 . OR . ND . GT . 3 65 ) THEN GOTO YEAR_OK 
$IF (ND . LE .MXDAY) THEN IRWRS=1 
$ IF ( IRWRS . EQ . 1 ) THEN GOTO OVERLAP 1 
$NO_OVERLAP 1 : 

$YYDDD= F$EXTRACT ( 0 , 2 , P2 ) + F$EXTRACT (0, 3, P3) 

$DIR/OUTPUT=DRECTRY . TMP ' DEFAULT_DIR' *' YYDDD' .DAT 
$IF $SEVERITY .EQ. 1 THEN GOTO GOOD_START 
$DELETE DRECTRY . TMP ; * 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "No files found with date ",' YYDDD' ,"! ” 

$ INQUIRE REDO "Want to reenter dates? [Y/N] " 

$IF REDO THEN GOTO NEW_YEAR 
$GOTO REMO VE_T AP E 
$ GOOD_S TART : 

SDELETE DRECTRY . TMP ; * 

$ISTART-ND 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT.- 

"If it is doubtful that all days will fit on this tape," 

$ WRITE SYS$OUTPUT "then try successive pairs of first and last days," 

$ WRITE SYS$OUTPUT "and put remaining days on a new tape." 

$ WRITE SYS$OUTPUT - 

"(It is OK to copy one day at a time, that is last day = first day.)" 
$YEAR_STILL_OK : 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE P3 "Last DAY to be copied? (in that same segment)" 
$LP3=F$LENGTH (P3) 

$IF (LP3 .GT . 3) THEN GOTO YE AR_S T I LL_OK 
$IF (LP3 .EQ . 2 ) THEN P3:= "0"P3'" 

$IF (LP3 . EQ . 1 ) THEN P3:= "00"P3'" 

$ND=F$INTEGER(F$EXTRACT(0,3,P3) ) 

$IF (NY.EQ. 80) .AND. (ND . LT . 1 . OR . ND . GT . 36 6 ) .OR. - 

(NY.NE.80) .AND. (ND . LT . 1 . OR . ND . GT . 3 65 ) THEN GOTO YEAR_STILL_OK 
$IF (ND . LT . I START) THEN GOTO WRONG_ORD E R_D AY S 
$IF (ND.GE.MNDAY) .AND. (IRWRS. EQ.l) THEN GOTO OVERLAP2 
$NO_OVERLAP 2 : 

$IRWRS=0 

$IF (ND.EQ. I START) THEN GOTO SINGLE_DAY 
$YYDDD= F$EXTRACT (0, 2, P2) + F$EXTRACT (0, 3, P3) 

$DIR/OUTPUT=DRECTRY . TMP ' DEFAULT_DIR' *' YYDDD' .DAT 
$IF SSEVERITY . EQ . 1 THEN GOTO GOOD_STOP 
SDELETE DRECTRY . TMP ; * 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYS$OUTPUT "No files found with date ",' YYDDD' ,"! " 

SINQUIRE REDO "Want to reenter both dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 

SINQUIRE REDO "Want to reenter last day? [Y/N]" 

$IF REDO THEN GOTO YE AR_S T I LL_OK 
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$GOTO REMO VE_T AP E 
$GOOD_STOP : 

$DELETE DRECTRY . TMP ; * 

$SINGLE_DAY: 

$ISTOP=ND 
$IDAY=ISTART 
$SET NOON 
$COPY_TO_TAPE : 

$IF (IDAY.LT.MNDAY) THEN MNDAY=IDAY 
$IF (IDAY.GT.MXDAY) THEN MXDAY=IDAY 

$ IF ( IDAY . LE . 9 ) THEN YYDDD=* F$STRING (NY) + "00" + F$STRING ( IDAY) 

$IF (IDAY.GT . 9) .AND. (IDAY. LE. 99) THEN - 
YYDDD= F$STRING (NY) + "0" + F$STRING ( IDAY) 

$ IF ( IDAY . GT .99) THEN YYDDD= F$STRING (NY) + F$STRING ( IDAY) 

$IF (IDAY. EQ. I START) THEN YYDDD 1 =YYDDD 
SCOPY/LOG 'DEFAULT_DIR' *' YYDDD' .DAT TAPE 
$IDAY=IDAY+1 

$IF (IDAY. LE . ISTOP) THEN GOTO COPY_TO_TAPE 
$NYOLD=NY 

$ WRITE SYSSOUTPUT " " 

$IF (ISTOP . EQ . I START) THEN WRITE SYSSOUTPUT - 
"Last segment copied: a single day ", 'YYDDD' 

$IF (ISTOP .GT . ISTART) THEN WRITE SYSSOUTPUT - 
"Last segment copied: days ",'YYDDD1'," through ", 'YYDDD' 

$ WRITE SYSSOUTPUT ” " 

$ INQUIRE MORE - 

"Are there more £ata files to be copied to this tape? [Y/N]" 

$IF MORE THEN GOTO NEW_YEAR 
$ WRITE SYSSOUTPUT " " 

$ INQUIRE GUIDE - 

"Do you want to add the Data User' s Guide at the end? [Y/N] " 

$IF GUIDE THEN - 

COPY/LOG 'DEFAULT_DIR' GUIDE. TAPE 
$TAPEDIR_LIS : 

$ WRITE SYSSOUTPUT " " 

SINQUIRE LISTING "Do you want a tape directory file (TAPEDIR.LIS) ? [Y/N]" 

$IF .NOT. LISTING THEN GOTO REMO VE_T AP E 

$DIR/SIZE/DATE/OUTPUT='DEFAULT_DIR' TAPEDIR.LIS TAPE 

SSET PROT= (G : RWED, W : RWED) TAPEDIR.LIS 

SGOTO REMO VE_T APE 

$WRONG_ORDER_DAYS : 

$ WRITE SYSSOUTPUT " " 

SWRITE SYSSOUTPUT "First and last days in wrong order!" 

SINQUIRE REDO "Want to reenter dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
SGOTO REMOVE_TAPE 
$WRONG_ORDER_YEARS : 

SWRITE SYSSOUTPUT " " 

SWRITE SYSSOUTPUT "Going backwards in years?" 

SINQUIRE SURE "Is this really the year you want? [Y/N]” 

$IF SURE THEN NYOLD=0 
$IF SURE THEN MNDAY=999 
$IF SURE THEN MXDAY=0 
$IF SURE THEN GOTO YEAR_OK 

SINQUIRE REDO "Want to correct year? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
SGOTO REMOVE TAPE 





$0VERLAP1 : 

$ WRITE SYSSOUTPUT ” " 

$ WRITE SYS$OUTPUT "Caution:" 

$ IF (ND . LT . MXDAY ) THEN WRITE SYS$OUTPUT - 
"First day precedes days already copied!" 

$IF (ND.EQ. MXDAY) THEN WRITE SYS$OUTPUT - 
"This day has already been copied!" 

S! IF (ND.EQ. MXDAY) THEN GOTO NEW_YEAR 
$ WRITE SYS$OUTPUT " " 

$ INQUIRE SURE "Are you sure there will be no duplication? [Y/N] 

$IF SURE THEN GOTO N0_0VERLAP1 

SINQUIRE REDO "Want to correct dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 
SGOTO REMO VE_T AP E 
SOVERLAP2 : 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT "Caution:" 

SIF(ND.GT.MNDAY) THEN WRITE SYSSOUTPUT - 
"Last day may be out of order!" 

$IF (ND.EQ. MNDAY) THEN WRITE SYSSOUTPUT - 
"Last day has already been copied!" 

$! IF (ND.EQ. MNDAY) THEN GOTO YEAR_STILL_OK 
$ WRITE SYSSOUTPUT " " 

SINQUIRE SURE "Are you sure there will be no duplication? [Y/N] 
$IF SURE THEN GOTO NO_OVERLAP2 

SINQUIRE REDO "Want to reenter both dates? [Y/N]" 

$IF REDO THEN GOTO NEW_YEAR 

SINQUIRE REDO "Want to reenter last day? [Y/N]" 

$IF REDO THEN GOTO YE AR_S T I LL_OK 
$ REMO VE_T APE : 

SSET ON 

SDISMOUNT 'DRIVE' 

SDEALLOCATE 'DRIVE' 

$ ! DELETE DRECTRY . TMP ; * 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "EXIT" 

$ WRITE SYSSOUTPUT "That's the end of this run." 

$ WRITE SYSSOUTPUT - 
"Tape is dismounted and 
SEXIT 


' ' FSLOGICAL ( "TAPE" ) ' is deallocated. 



$ | ************** LOCKHEED PLASMA COMPOSITION EXPERIMENT ON ISEE-1 ************** 
$ ! ****** DELETE_ISEEMS.COM ****** 

$ ! 

$! Interactive command file for deleting ISEEMS archival files from disk 
$! (can transfer the deleting to a batch job via REMOVE_ISEEMS__FILES . COM) 

$ ! 

$! by 0. W. Lennartsson 

$! Lockheed Martin Missiles & Space 

$ | Advanced Technology Center 

$! Dept. Hl-11, Bldg. 255 

$! 3251 Hanover Street 

$! Palo Alto, CA 94304 

$ ! 

$! Telephone (650) 424-3259 
$ ! 

$! Developed under NASA Contract NAS5-33047 and NASA Grant NAGW-4496 
$ ! 

$! Operating system used: VAX/VMS version V5 . 5 
$ ! 

$ ! To execute this file, simply enter @DELETE_ISEEMS, if the file resides in the 
$! same directory used for the data files, otherwise specify directory name in 
$! the command: @directory_nameDELETE_ISEEMS (your default directory should be 
$! the one used for data files) . 

<; I ****************************************************************************** 

$ WRITE SYS$OUTPUT " " 

$DEFAULT_D I R = F$LOGICAL ( "SYSSDISK" ) +F$DIRECTORY ( ) 

SIThis directory may or may not be used for the command files. 

$!If not, change next command line to appropriate directory: 

$COMMAND_D I R = DEFAULT_DIR 
$ ! 

$ WRITE SYS$OUTPUT "Your default directory is: ", " DEFAULT_DIR' 

$ WRITE SYS$OUTPUT "Is this where the ISEEMS data files are?" 

$ WRITE SYSSOUTPUT "If not, Ctrl/Y and reset." 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "Files are deleted by date, starting with the 
"first day entered, " 

$ WRITE SYSSOUTPUT "and ending with the second day 
"-assuming one and the same year." 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYSSOUTPUT "If files run into second year, ",- 
"repeat this operation for that year's days." 

$ WRITE SYSSOUTPUT " " 

SINQUIRE P2 "What YEAR? (two digits)” 

$LENGTH_OF_YEAR : 

$LP2=F$LENGTH (P2) 

$IF (LP2 . EQ . 2 ) THEN GOTO CHECK_YEAR 

SINQUIRE P2 "Enter two digits for the year: 77 through 82" 

SGOTO LENGTH_OF_YEAR 
$WRONG_YEAR : 

SINQUIRE P2 "Enter either of 77, 78, 79, 80, 81, or 82" 

$CHECK_YEAR : 

$NY=F$ INTEGER (F$EXTRACT (0, 2, P2) ) 

$IF (NY . LT .77. OR . NY . GT . 82 ) THEN GOTO WRONG_YEAR 
$YEAR_OK : 

$ WRITE SYSSOUTPUT " " 

SINQUIRE P3 "First DAY? (one to three digits)" 


/O 



$LP3=F$ LENGTH (P3) 

$IF (LP3 .GT . 3 ) THEN GOTO YEAR_OK 
$IF (LP3 . EQ . 2.) THEN P3:= ”0"P3'" 

$IF(LP3.EQ.l) THEN P3:= "00"P3'” 

$ND=F$INTEGER (F$EXTRACT ( 0 , 3,P3) ) 

$IF (NY . EQ .80) .AND. (ND . LT . 1 . OR . ND . GT . 3 6 6 ) .OR. - 

(NY.NE.80) .AND. (ND . LT . 1 . OR . ND . GT . 3 65 ) THEN GOTO YEAR_OK 
$YYDDD= F$EXTRACT ( 0 , 2 , P2 ) + F$EXTRACT (0, 3, P3) 

$DIR/OUTPUT=DRECTRY . TMP ' DEFAULT_DIR' *' YYDDD' .DAT 
$IF $SEVERITY . EQ . 0 THEN GOTO NO_FILE 
SDELETE DRECTRY . TMP ; * 

$ISTART=ND 
$YEAR_STILL_OK : 

$ WRITE SYS$OUTPUT ” " 

SINQUIRE P3 "Second DAY? (last day of same year to be deleted)" 
$LP3=F$LENGTH (P3) 

$IF (LP3 .GT . 3) THEN GOTO YEAR_STILL_OK 
$IF (LP3 . EQ . 2 ) THEN P3:= "0"P3'" 

$IF (LP3 . EQ . 1) THEN P3:= ”00"P3'" 

$ND=F$ INTEGER (F$EXTRACT (0, 3,P3) ) 

$IF (NY.EQ. 80) .AND. (ND . LT . 1 . OR . ND . GT . 3 66 ) .OR. - 

(NY.NE.80) .AND. (ND . LT . 1 . OR . ND . GT . 3 65 ) THEN GOTO YE AR_S TIL L_0 K 
$YYDDD= F$EXTRACT (0, 2, P2) + F$EXTRACT (0, 3, P3) 

$DIR/OUTPUT=DRECTRY . TMP ' DEFAULT_DIR' *' YYDDD' .DAT 
$IF $SEVERITY . EQ . 0 THEN GOTO NO_FILE 
$DELETE DRECTRY . TMP ; * 

$ISTOP=ND 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "If you are deleting a large number of files it may ",- 
"take a few minutes . " 

$WRITE SYS$OUTPUT "In the meantime this terminal will be tied up, 

"unless you transfer this job" 

$ WRITE SYSSOUTPUT "to a batch queue." 

$ WRITE SYS$OUTPUT " " 

$ INQUIRE BAT "Would you rather delete in BATCH? [Y/N] " 

$IF .NOT. BAT THEN GOTO NO_BAT 

SINQUIRE BATQU "Enter name of batch queue" 

$SUBMIT/QUEUE=' BATQU' /LOG_FILE=' DEFAULT_DIR' /NOPRINT- 
/PARAM=('DEFAULT_DIR, 'NY, ' ISTART, ' ISTOP) - 
' COMMAND_DIR' REMOVE_ISEEMS_FILES 
$ WRITE SYS$OUTPUT " " 

$ WRITE SYSSOUTPUT "This terminal is yours again" 

$ WRITE SYS$OUTPUT ” ” 

$EXIT 
$NO_BAT : 

$ WRITE SYSSOUTPUT " " 

$ WRITE SYS$OUTPUT "As you wish!" 

$ WRITE SYSSOUTPUT " " 

$IDAY=ISTART 
$DATA_DELETE : 

$IF (IDAY.LE. 9) THEN YYDDD= F$STRING (NY) + "00" + F $ STRING ( IDAY) 

$IF ( IDAY . GT . 9) .AND. (IDAY.LE. 99) THEN - 
YYDDD = FSSTRING (NY) + "0" + F$STRING ( IDAY) 

$ IF ( IDAY . GT .99) THEN YYDDD= FSSTRING (NY) + F $ STRING ( IDAY) 

SDELETE ' DEFAULT_DIR' *' YYDDD' .DAT; * 

$IDAY=IDAY+1 

$IF ( IDAY. GT. ISTOP) THEN GOTO FILES_LEFT 


// 



$GOTO DATA_DELETE 
$FILES_LEFT: 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "Files left:” 

$ WRITE SYSSOUTPUT ” " 

$DIR/SIZE=ALL/DATE/PROT ED* , EX* , MC* , MD* , MS * 

$EXIT 

$NO_FILE: 

$ WRITE SYS$OUTPUT " " 

$ WRITE SYS$OUTPUT "Files missing: No ED", ' YYDDD' , " .DAT, 
$ WRITE SYS$OUTPUT "Check dates and reenter command." 
SDELETE DRECTRY . TMP ; * 

$EXIT 


/X~ 


etc . 



$;************** LOCKHEED PLASMA COMPOSITION EXPERIMENT ON ISEE-1 ************** 
$ ! ****** REMOVE_ISEEMS_FILES.COM ****** 

$ ! 

$! Batch command file for deleting ISEEMS archival files from disk 

$! (must be submitted from the interactive DELETE_ISEEMS.COM) 

$ ! 

$! by 0. W. Lennartsson 

$ ! Lockheed Martin Missiles & Space 

$ \ Advanced Technology Center 

$ ! Dept. Hl-11, Bldg. 255 

$! 3251 Hanover Street 

$! Palo Alto , CA 94304 

$ ! 

$! Telephone (650) 424-3259 

$! 

$! Developed under NASA Contract NAS5-33047 and NASA Grant NAGW-4496 

$ ! 

$! Operating system used: VAX/VMS version V5 . 5 

c;l****************************************************************************** 

$DEFAULT_DIR=P1 

$SET DEFAULT ' DEFAULT_DIR' 

$NY=F $ INTEGER ( P 2 ) 

$ ISTART=F$ INTEGER ( P 3 ) 

$ISTOP=F$INTEGER(P4) 

$SET NOON 
$SET NOVERIFY 
$IDAY=ISTART 
$DATA_DELETE : 

$ IF ( IDAY . LE . 9 ) THEN YYDDD= F$STRING (NY) + "00" + F$STRING ( IDAY) 

$IF (IDAY.GT. 9) .AND. (IDAY. LE. 99) THEN - 
YYDDD= F$STRING (NY) + "0" + F$STRING ( IDAY) 

$IF (IDAY.GT. 99) THEN YYDDD= F$STRING (NY) + F$STRING ( IDAY) 

$DELETE * ' YYDDD ' . DAT ; * 

$IDAY=IDAY+1 

$IF ( IDAY . LE . ISTOP ) THEN GOTO DATA_DELETE 

$SET VERIFY 

$EXIT 
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Appendix 5 


Sample Magnetic Tape Directory Listing 

(1 page) 
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